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Abstract. In this article, we derive a Stratonovich and Skorohod type change of vari- 
ables formula for a multidimensional Gaussian process with low Holder regularity 7 (typ- 
ically 7 < 1/4). To this aim, we combine tools from rough paths theory and stochastic 
analysis. 



1. Introduction 

Starting from the seminal paper [7], the stochastic calculus for Gaussian processes has 
been thoroughly studied during the last decade, fractional Brownian motion being the 
main example of application of the general results. The literature on the topic includes 
the case of Volterra processes corresponding to a fBm with Hurst parameter H > 1/4 (see 
[TJ H2]), with some extensions to the whole range H G (0, 1) as in |2J [HI ITTj . It should be 
noticed that all those contributions concern the case of real valued processes, this feature 
being an important aspect of the computations. 

In a parallel and somewhat different way, the rough path analysis opens the possibility 
of a pathwise type stochastic calculus for general (including Gaussian) stochastic pro- 
cesses. Let us recall that this theory, initiated by T. Lyons in [20] (see also |HJ ETJ [T3] for 
introductions to the topic), states that if a 7-Holder process x allows to define sufficient 
number of iterated integrals then: 

(1) One gets a Stratonovich type change of variable for f(x) when / is smooth enough. 

(2) Differential equations driven by x can be reasonably defined and solved. 

In particular, the rough path method is still the only way to solve differential equations 
driven by Gaussian processes with Holder regularity exponent less than 1/2, except for 
some very particular (e.g. Brownian, linear or one-dimensional) situations. 

More specifically, the rough path theory relies on the following set of assumptions: 

Hypothesis 1.1. Let 7 G (0,1) and x : [0,T] R d be a ^-Holder process. Consider 
also the n th order simplex S ny T = {{u\, ■ ■ ■ , u n ) : < u\ < ■ ■ • < u n < T} on [0, T\. 
The process x is supposed to generate a rough path, which can be understood as a stack 
{x 11 ; n < |_1/7_|} of functions of two variables satisfying the following three properties: 

(1) Regularity: Each component o/x n is wy-Holder continuous (in the sense of the Holder 
norm introduced in ^^)) for all n < [1/tJ; an d x ^ = x t — x s . 
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(2) Multiplicativity: Letting (5x. n ) sut := — x" u — x" f for (s,u,t) £ Sz t T, one requires 

n-l 

(dX- n )sut(h, ■ ■ ■ ,i n ) — 2^ X "w( Zl ' • • • '^m) X «t ni (^m+l) • • • ) z 'n)- (1) 

ni=l 

Geometricity: For any n, m such that n + m < [1/tJ ( s ; £ <^2,T, we /icroe: 

X si(^l' • • • 5 i"n) x s t (jl) • • • ) jm) = x s t + (&1 , • • • > k n+m ) , (2) 

keSh(i,j) 

where, for two tuples t,J, stands for the set of permutations of the indices contained 

in and Sh(i,j) is a subset o/E^j) defined by: 

Sh(i,j) = {cr £ a does not change the orderings oft and]} . 

With this set of abstract assumptions in hand, one can define integrals like J f(x) dx 
in a natural way (as recalled later in the article), and more generally set up the basis of a 
differential calculus with respect to x. Notice that according to T. Lyons terminology |21| . 
the family {x 11 ; n < [1/ -yj } is said to be a weakly geometric rough path above x. 

Without any surprise, some substantial efforts have been made in the last past years 
in order to construct rough paths above a wide class of Gaussian processes, among which 
emerges the case of fractional Brownian motion. Let us recall that a fractional Brownian 
motion B with Hurst parameter H £ (0,1), defined on a complete probability space 
(Jl, J 7 , P), is a d-dimensional centered Gaussian process. Its law is thus characterized by 
its covariance function, which is given by 

E [B t {i)B s {i)\ = l - (t 2H + s 2H -\t- s\ 2H ) l (i=i) , 8 ,te R+. (3) 

The variance of the increments of B is then given by 

E [(B t (i) - B s (t)) 2 ] =(t- s) 2H , (s, t) £ S 2 , T , i = l,...,d, 

and this implies that almost surely the trajectories of the ffim are 7-Holder continuous for 
any 7 < H. Furthermore, for H = 1/2, fBm coincides with the usual Brownian motion, 
converting the family {B = B H ; H £ (0, 1)} into the most natural generalization of this 
classical process. This is why B can be considered as one of the canonical examples of 
application of the abstract rough path theory. 

Until very recently, the rough path constructions for fBm were based on pathwise type 
approximations of B, as in (31 [231 EH]. Namely, these references all use an approximation 
of B by a regularization B e , consider the associated (Riemann) iterated integrals B n,e and 
show their convergence, yielding the existence of a geometric rough path above B. These 
approximations all fail for H < 1/4. Indeed, the oscillations of B are then too heavy 
to define even B 2 following this kind of argument, as illustrated by [S]. Nevertheless, 
the article [22] asserts that a rough path exists above any 7-Holder function, and the 
recent progresses |2B| [25] show that different concrete rough paths above fBm (and more 
general processes) can be exhibited, even if those rough paths do not correspond to a 
regularization of the process at stake. 

Summarizing what has been said up to now, there are (at least) two ways to handle 
stochastic calculus for Gaussian processes: (i) Stochastic analysis tools, mainly leading 
to a Skorohod type integral (ii) Rough paths analysis, based on the pathwise convergence 
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of some Riemann sums and giving rise to a Stratonovich type integral. Though some 
efforts have been made in [3] in order to relate the two approaches (essentially for a fBm 
with Hurst parameter H > 1/4), the current article proposes to delve deeper into this 
direction. Namely, we plan to tackle three different problems: 

(1) We show that, starting from a given rough path of order N above a <i-dimensional 
process x, one can derive a Stratonovich change of variables of the form 

f(x t ) - f(x s ) = ^2 9 ifM dx u(i) ■= Jst (V/(ar u ) dx u ) , (4) 
i=i Js 

for any / £ C N+1 (M. d \ M), and where dif stands for df/dxi. This formula is not new, 
and is in fact an immediate consequence of the powerful stability theorems which can be 
derived from the abstract rough paths theory (see e.g [9]). However, we have included 
these considerations here for several reasons: (i) This paper not being dedicated to rough 
paths specialists, we find it useful to include a self contained, short and simple enough 
introduction to equation (J4]) (ii) Our proof is slightly different from the original one, in 
the sense that we only rely on the algebraic and analytic assumptions of Hypothesis 11.11 
rather than on a limiting procedure (Hi) Proving (J4J) is also a way for us to introduce all 
the objects and structures needed later on for the Skorohod type calculus. In particular, 
we derive the following representation for the integral i7st(V/(x u ) dx u ): consider a family 
of partitions U st = {s = t , . . . , t n = t} of [s,t], whose mesh tends to 0. Then, denoting 

byiv= 

n-l N-l , 

J st (Vf(x u ) dx u ) = lim £ £ yd^JM *} qtq+1 (i k ) ■ ■ ■ x\ tq+1 {h) (5) 

3t ' q=0 k=0 

These modified Riemann sums will also be essential in the analysis of Skorohod type 
integrals. 

(2) We then specialize our considerations to a Gaussian setting, and use Malliavin calculus 
tools (in particular some elaborations of E])- Namely, supposing that x is a Gaussian 
process, plus mild additional assumptions on its covariance function, we are able to prove 
the following assertions: 

(i) Consider a C 2 (M. d ; M.) function / with exponential growth, and < s < t < oo. Then 
the function u t— > lr S) t)(u)V/(a; w ) lies into the domain of an extension of the divergence 
operator (in the Malliavin calculus sense) called 5°. 

(ii) The following Skorohod type formula holds true: 

f(x t ) - f(x s ) = 5" (1 M V/(s)) + \ J Af(x u ) R' u du, (6) 

where A stands for the Laplace operator, u t— > R u := E[|x u (l)| 2 ] is assumed to be a 
differentiable function, and R' stands for its derivative. 

It should be emphasized here that formula is obtained by means of stochastic analysis 
methods only, independently of the Holder regularity of x. Otherwise stated, as in many 
instances of Gaussian analysis, pathwise regularity can be replaced by a regularity on 
the underlying Wiener space. When both, regularity of the paths and on the underlying 
Wiener space, are satisfied we obtain the relation between the Stratonovich type integral 
and the extended divergence operator. 
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Let us mention at this point the recent work [TH] that considers similar problems as 
ours. In that article, the authors define also an extended divergence type operator for 
Gaussian processes (in the one-dimensional case only) with very irregular covariance and 
study its relation with a Stratonovich type integral. For the definition of the extended 
divergence, some conditions on the distributional derivatives of the covariance function R 
are imposed, one of them being that d1 t R st satisfies that Jx{ds, dt) := dg t R st (t — s) (that is 
well defined) is the difference of two Radon measures. Our conditions on R are of different 
nature, we suppose more regularity but only for the first partial derivative of R and the 
variance function. On the other hand, the definition of the Stratonovich type integral 
in [19] is obtained through a regularization approach instead of rough paths theory. As 
a consequence, some additional regularity conditions on the Gaussian process have to be 
imposed, while we just rely on the existence of a rough path above x. 

(3) Finally, one can relate the two stochastic integrals introduced so far by means of 
modified Wick-Riemann sums. Indeed, we shall show that the integral 5° (l[s,t) 
introduced at relation ([6]) can also be expressed as 

n-l N-l , 

? = |n iim o E E y€Xif(^)^ qtq+1 (i k )<>- ■ ■ ^(iOox^co, (7) 

st| g=0 fc=0 

where the (almost sure) limit is still taken along a family of partitions U st = {s = 
to, . . . ,t n = t} of [s,t] whose mesh tends to 0, and where o stands for the usual Wick 
product of Gaussian analysis. This result can be seen as the main contribution of our 
paper, and is obtained by a combination of rough paths and stochastic analysis methods. 
Specifically, we have mentioned that the modified Riemann sums in ([5]) can be proved 
to be convergent by means of rough paths analysis. Our main additional technical task 
will thus consist in computing the correction terms between those Riemann sums and the 
Wick-Riemann sums which appear in (J7J). This is the aim of the general Proposition 16.71 
on Wick products, which has an interest in its own right, and is the key ingredient of our 
proof. It is worth mentioning at this point that Wick products are usually introduced 
within the landmark of white noise analysis. We rather rely here on the introduction 
given in [T7], using the framework of Gaussian spaces. Let us also mention that Riemann- 
Wick sums have been used in [ L 8J to study Skorohod stochastic calculus with respect to 
(one-dimensional) fBm for H greater than 1/2, the case of 1/4 < H < 1/2 being treated 
in [27] . We go beyond these case in Theorem 16. 8[ and will go back to the link between 
our formulas and the one produced in [27] at Section 16.31 

In conclusion, this article is devoted to show that Stratonovich and Skorohod stochas- 
tic calculus are possible for a wide range of Gaussian processes. A link between the 
integrals corresponding to those stochastic calculus is made through the introduction of 
Riemann- Wick modified sums. On the other hand, the reader might have noticed that 
the integrands considered in our stochastic integrals are restricted to processes of the form 
V/(x). The symmetries of this kind of integrand simplify the analysis of the Stratonovich- 
Skorohod corrections, reducing all the calculations to corrections involving x 1 only. An 
extension to more general integrands would obviously require a lot more in terms of Wick 
type computations, especially for the terms involving x k for k > 2, and is deferred to a 
subsequent publication. 
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Here is how our paper is organized: Section [2] recalls some basic elements of rough paths 
theory which will be useful in the sequel. Then, as a warmup for the non initiated reader, 
we derive a Stratonovich change of variable formula in the case of a rough path of order 
2 at Section |3j The case of a rough path of arbitrary order is then treated at Section HI 
We obtain a Skorohod change of variable with Malliavin calculus tools only at Section |5j 
Finally, the representation of this Skorohod integral by Wick-Riemann sums is performed 
at Section 

2. Some elements of algebraic integration 

As already mentioned in the introduction, our stochastic calculus will appeal to the 
algebraic integration theory, which is a variant of the rough paths theory introduced 
in [13], and for which we also refer to |15| for a detailed introduction. 

2.1. Increments. The extended pathwise integration we will deal with is based on the 
notion of 'increments', together with an elementary operator 5 acting on them. The 
algebraic structure they generate is described in [131 fT5] . but here we present directly the 
definitions of interest for us, for sake of conciseness. First of all, for an arbitrary real 
number T > 0, a vector space V and an integer k > 1 we denote by Ck{V) the set of 
functions g : [0, T] k — > V such that Qt^-t k — whenever tj = for some % < k — 1. Such 
a function will be called a (k — l)-increment, and we set C*{V) = U k >iC k (V). We can now 
define the announced elementary operator 5 on C k (V): 

fe+i 

5 : C k (V) C k+1 (V), (6g) tl ... tk+1 = £(-1) fc -V/ tl ... {< ... tfr ., : (8) 

i=l 

where ti means that this particular argument is omitted. A fundamental property of 5, 
which is easily verified, is that 55 = 0, where 55 is considered as an operator from Ck{V) 
to C k+2 {V). We denote ZC k {V) = C k (V) D Ker5 and BC k (V) = C k (V) n Im5. 

Some simple examples of actions of 5, which will be the ones we will really use through- 
out the paper, are obtained by letting g G C\ and h G Then, for any s,u,t G [0, T], 
we have 

(8g)st = 9t- 9s, and {5h) sut = h st - h su - h ut . (9) 

Furthermore, it is easily checked that ZC k+ i(V) = BC k (V) for any k > 1. In particular, 
the following basic property holds: 

Lemma 2.1. Let k > 1 and h G ZC k+ i(V). Then there exists a (non unique) f G C k (V) 
such that h = 5f . 

Proof. This elementary proof is included in [T3j, and will be omitted here. However, let 
us mention that /^...^ = (— 1) 1 hot 1 ...t h is a possible choice. □ 

Observe that Lemma 12.11 implies that all the elements h G C2{V) such that 5h = 
can be written as h = 5f for some (non unique) / G Ci(V). Thus we get a heuristic 
interpretation of 5\c 2 (v)'- h measures how much a given 1-increment is far from being an 
exact increment of a function, i.e., a finite difference. 
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Notice that our future discussions will mainly rely on /c-increments with k < 2, for 
which we will make some analytical assumptions. Namely, we measure the size of these 
increments by Holder norms defined in the following way: for / G C^iy) let 

ix = sup J^L an d (%(V) = {/ G C 2 (V); < 00} . (10) 



s,te[o,T] 



t - s 



Obviously, the usual Holder spaces (V) will be determined in the following way: for a 
continuous function g G Ci(V), we simply set 

NUHNU (ii) 

and we will say that g G Cf (V) iff ||p || ^ is finite. Notice that || • || M is only a semi-norm on 
Ci(V), but we will generally work on spaces of the type 

Ci, a (V) = {g : [0,71 ^V;g = a, \\g\\, < 00}, (12) 

for a given a G V, on which || (7 1| M defines a distance in the usual way. For h G Cs(V) set 
in the same way 



,., :t ...n (1 3 ) 



7,/o 



s,u,tG[0,T] 



kt - spit - u\p 



where the last infimum is taken over all sequences {hi G Cs(V)} such that h = hi and 
for all choices of the numbers G (0, z). Then || • is easily seen to be a norm on C 3 (V), 
and we set 

C%(V) :={heC 3 (V); ||%<oo}. 

Eventually, let C\ + (V) = U^iC^V), and notice that the same kind of norms can be 
considered on the spaces ZCs(V), leading to the definition of some spaces ZC^iV) and 
ZC\ + {V). 

With these notations in mind the following proposition is a basic result, which belongs 
to the core of our approach to pathwise integration. Its proof may be found in a simple 
form in |15] . 

Proposition 2.2 (The A-map). There exists a unique linear map A : ZC\ + {V) — > C\ + {V) 
such that 

6 A = Id zc i+( V j and A5 = Id c i+^ v y 

In other words, for any h G C\ + {V) such that 5h = there exists a unique g = A(h) G 
C] + (V) such that Sg = h. Furth ermore, for any ji > 1, the map A is continuous from 
ZC^iV) to C% (V) and we have 

||A%<^_||%, h G ZC^iV). (14) 

Let us mention at this point a first link between the structures we have introduced so 
far and the problem of integration of irregular functions. 
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Corollary 2.3. For any 1-increment g G C<z(V) such that Sg G C\ + , set Sf = (Id — AS)g. 
Then 

n-l 

( 5 f)st = ,J i V 1 n Xl^**+i' 

list] — ^o . 

where the limit is over any partition U st = {to = s, . . . ,t n = t} of [s, t], whose mesh tends 
to zero. Thus, the 1-increment Sf is the indefinite integral of the 1-increment g. 

Proof. Just consider the equation g = Sf + ASg and write 

n—l n—X n— 1 

Sn st = ^Quu+i = ^2($f)tit i+1 + ^2(ASg) uu+1 

i=0 i=0 i=0 

n-l 

= (*/)*+ x>fow 

i=0 

Then observe that, due to the fact that A5<? G Cg 4 "^), the last sum converges to zero. 

□ 

2.2. Computations in C*. Let us specialize now to the case V = R, and just write 
for (R). Then (C*, 5) can be endowed with the following product: for g G C n and h G C m 
let gh be the element of C n+m _i defined by 

(gh) tl ,...,t m+n+1 = 9t u ...,t n ht n ,...,t m+n - 1 , h, ■ ■ ■ ,tm+n-i G [0,T]. (15) 
In this context, we have the following useful properties. 

Proposition 2.4. The following differentiation rules hold true: 

(1) Let g G C\ and h G C\. Then gh G C\ and 

S(gh) = Sgh + gSh. (16) 

(2) Let g G Ci and h £ C 2 - Then gh G C2 and 

S(gh) = Sgh-gSh. (17) 
Lei g G C2 and h E C±. Then gh G C2 and 

S(gh) = Sgh + gSh. (18) 
Proof. We will just prove (1161) . the other relations being just as simple. If g,h G C\, then 

i 5 (9h)] st = g t ht - g s h s = g s (h t - h s ) + (g t - g s ) h t = g s (Sh) st + (Sg) st h t , 
which proves our claim. 

□ 

The iterated integrals of smooth functions on [0, T] are particular cases of elements of 
C%, which will be of interest 

for us. Let us recall some basic rules for these objects: consider / G and g G Cf 3 , 
where Cf 3 denotes the set of smooth functions on [0, T\. Then the integral j f dg, which 
will be denoted indistinctly by J f dg or J(f dg), can be considered as an element of 
Namely, for (s,t) G <S2,r we set 



JMdg) = 
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The multiple integrals can also be defined in the following way: given a smooth element 
h G C|° and (s,t) G <5>2,t, we set 



J s t(hdg) 




In particular, for f l G Cf, / 2 G Cf 3 and / 3 G Cf the double integral Jst(Pdpdp) is 
defined as ^ 

Now suppose that the nth order iterated integral of f n+1 df n ■ ■ ■ df 2 , which is denoted by 
J(f n+1 df n ■ ■ ■ df 2 ), has been defined for p G Cf. Then, if f 1 G Cf, we set 

JM n+1 df n ■ ■ ■ dpdf 1 ) = f J su {f n+l dp ■ ■ ■ dp) dp, (19) 

J s 

which recursively defines the iterated integrals of smooth functions. Observe that an nth 
order integral J{df n ■ ■ ■ dpdf 1 ) can be defined along the same lines, starting with 

J(df) = Sf, 

J st (df 2 dp)= [ J su (df 2 )dp u = [ (sp) su dp, 

J s J s 

and so on. 

The following relations between multiple integrals and the operator 5 will also be useful. 
The reader is sent to [12] for its elementary proof. 

Proposition 2.5. Let f G Cf and g G Cf 3 . Then it holds that 

5g = J(dg), 5 (J(fdg)) = 0, 5 (J(dfdg)) = (5f)(5g) = J(df)J(dg), 

and 

n-l 

5 {J {df n ■ ■ ■ df 1 )) = J {df n ■ ■ ■ dp +1 ) J {dp ■ ■ ■ df 1 ) . 

i=l 

3. STRATONOVICH CALCULUS OF ORDER 2 

This section is devoted to establish an Ito-Stratonovich change of variable formula for 
a process x G Cj(M. d ), with 1/3 < 7 < 1/2, provided this process generates a (weakly 
geometric) Levy area. It is intended as a warm up for the general change of variable of 
the next section, especially for those readers who might not be acquainted to rough paths 
techniques. 

3.1. Weakly controlled processes. Recall that we have in mind to give a change of 
variable formula for f(x) when a; is a function in Cj(R d ) with 7 > 1/3 and / is a sufficiently 
smooth function. In this case, the rough path above x is reduced to a second order iterated 
integral, and the multiplicative property (JTJ of the path can be read as: 

Hypothesis 3.1. The path x is W 1 -valued ^-Holder with 7 > 1/3 and admits a Levy area, 
that is a process x 2 G C^iW 14 ) satisfying 

5x 2 = x 1 <g> x 1 , %. e. [(Sx 2 ) sut \ (i,j) = [x 1 ^)]^ 1 ^)]^, 

for s,u,t G Ss t T and i,j G {1, . . . , d}. 
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Let us now be more specific about the global strategy we will adopt in order to obtain 
our Stratonovich type formula. First of all, we shall define integrals with respect to x for 
a class of integrands called weakly controlled processes, that we proceed to define. Notice 
that in the following definition we use for the first time the convention of summation over 
repeated indices, which will prevail until the end of Section HI 

Definition 3.2. Let z be a process in Cj(R n ) with 1/3 < 7 < 1/2 (that is, N : = [1/tJ = 
2). We say that z is a weakly controlled path based on x and starting from a if zq = a, 
which is a given initial condition in R n , and 5z £ C\ (R n ) can be decomposed into 

Sz(i) = C(Mi) xl («i) + r{i), i-e. {6z{i)) st = «i) x st(*i) + r st{i), (20) 
for all (s,t) £ S2,t- In the previous formula, we assume £ £ Cj(R n ' d ), and r is a regular 
part such that r £ C^iR 71 ). The space of weakly controlled paths starting from a will be 
denoted by Q lta (R n ), and a process z £ Q 7ia (M n ) can be considered in fact as a couple 
(z,Q- The natural semi-norm on Q K>a (M. k ) is given by 

M[z- Q 7 , a (M n )] = AT[z; Cj(R n )} + A/1C; C™(M. n ' d )] + Cj{R n ' d )} + A/>; C 2 2 \R r % 
withAf[g;C?] defined by fOP and (V)) = sup < s < T \( s \ v . 

With this definition at hand, we will try to obtain our change of variables formula in 
the following way: 

(1) Study the decomposition of f(x) as weakly controlled process, when / is a smooth 
function. 

(2) Define rigorously the integral J z u dx u = J{zdx) for a weakly controlled path z 
and compute its decomposition f l20|) . 

(3) Compare the decompositions of f(x) and J V/(x) dx, and show that they coincide, 
up to a term with Holder regularity greater than 1. 

In this section, we will concentrate on the first point of the program. 

Let us see then how to decompose f(x) as a controlled process when / is a smooth 
enough function, for which we first introduce a convention which will hold true until the 
end of the paper: for any smooth function / : R d — > R, k > 1, (ji, . . . , if.) £ {1, . . . , d} k 
and £ £ R d , we set 

With this notation in hand, our decomposition result is the following: 

Proposition 3.3. Let f : R d — > R be a C% function such that f(xo) = a, and set z = f(x) . 
Then z £ Q 7 , a ; an d it can be decomposed into 5z = (5x + r, with 

C{i)=dif{x) and r = 5f(x) - dif(x) x 1 (z). 

Furthermore, 

M[z- Q 7 , Q ] < c ft T (1 + Af 2 [x;Cj(R d )}) ■ (22) 
Proof. The algebraic part of the assertion is straightforward. Just write 

{SzU = f(x t ) - f(x s ) = dj(x s ) x^z) + r st , (23) 
which is the desired decomposition. 
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In order to give an estimate for J\f[z; Q 7)0 (M n )], one has of course to establish bounds 
for Af[z;Cj{R n )), ^f[C; CJ(R% Af[(; C^°(R d )] and A/>;C 2 27 ]. These estimates are readily 
obtained from decomposition f )23|) . and details are left to the reader. 

□ 

Remark 3.4. The algebraic part of the above proposition remains true if we only suppose 
that / G C 2 (R d ). Indeed, since / together with its first and second order partial derivatives 
and x are continuous functions on a compact set, we have that ((i) = dif(x) G C\ and 
r = 8f(x) — dif(x)~x}{i) G C\. Nevertheless, the inequality norm (1221) fails and M[z\ Q 7)0 ] 
cannot be bounded by terms only depending on the Holder norm of x. 

3.2. Integration of weakly controlled paths. Let us now turn to the integration of 
weakly controlled paths, which is summarized in the following theorem. 



Theorem 3.5. For a given 1/3 < 7 < 1/2, let x be a process satisfying Hypothesis \3.1[ 
Furthermore, let m G Q ly b(R d ) with decomposition itlq = b G R d and 

5m{%) = /x(i,ii) x^ii) + r(i), where // G Cj(R d ' d ), r G C^{R d ). (24) 
Define z by zq = a G R and 

Sz = m(i) x 1 (z) + ii) x 2 (z 1; i) — A (r(i) x 1 (i) + 5[i(i, z'i) x 2 (z 1; z)) . (25) 
Finally, set 

J st (mdx) = J (m u , dx u )^ d = (Sz) st . 

Then: 

(1) z is well-defined as an element of Q 7)0 (R) ; and coincides with the Riemann-Stieltjes 
integral of z with respect to x whenever these two functions are smooth. 

(2) The semi-norm of z in Q 7|0 (M.) can be estimated as 

M[z\ Q 7 ,a(K)] < c x (1 + M[m- Q 7 , fe (M d )]) , (26) 
for a positive constant c x which can be bounded as follows: 

c x <c (A/lx 1 ; C](R d )} + W[x 2 ; CpiR^)]) , for a universal constant c. 

(3) It holds 



n-l 

Jstimdx) = lim V \m tq {i) xj, )t , +1 (i) + z '0 ^t^u i) , (27) 

q=0 



for any < s < t < T, where the limit is taken over all partitions U st = {s = t , . . . , t n 
t} of[s,t], as the mesh of the partition goes to zero. 



Before going into the technical details of the proof, let us see how to recover (125)) in the 
smooth case, in order to justify our definition of the integral. (Notice however that (l2Tj) 
corresponds to the usual definition in the rough paths theory [21], which gives another 
kind of justification.) 

Let us assume for the moment that x is a smooth function and that m G CJ°(K. d ) admits 
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the decomposition ((Ml) with \i G C%°(R d > d ) and r G C|°(R d ). Then f*{m u , dx u ) is well- 
defined, and we have 

t i-t 
(m u: dx u } Rd = m s (i) [x t (i) - x s (i)] + / [m u {i) - m s (i)] dx u (i) 

J s 

for s < t, respectively 

J{mdx) = m(i) x 1 (i) + J (Sm(i) dx(i)) . 
Let us now plug the decomposition f[24]) into this expression, which yields 

J st (mdx) = m s (i) x^(z) + / /i s (i, i x ) x^^) dx u (i) + J(r dx) 

J s 

= m s (i)x] t (i) +/x s (z,zi)Xg t (zi,z) + J(rdx). (28) 

Notice that the terms mSx and /ix 2 in (128]) are well-defined as soon as x and x 2 are 
defined themselves. In order to push forward our analysis to the rough case, it remains 
to handle the term J{r dx). Thanks to (128]) we can write 

J{r dx) = J{m dx) — m(i) x 1 ^) — /i(z, i±) x 2 («i, i), 

and let us analyze this relation by applying 5 to both sides. Using the second part of 
Proposition 12.41 and the Proposition 12.51 yields 

5[J(rdx)} = — 5 [m(i) x 1 (i)] — S [fi(i, i\) x 2 (zi, ij\ 

= —Sm(i) x 1 (i) — 8fi(i, ii) x 2 (zi, i) + /i(z, i\) x 1 (ii) x 1 (i) 

= - [/itMijx 1 ^) +r{i)] x*(i) - Sfifaijypfai) + //(i, z'i) x^z'i) x 1 ^) 

= -^(m'Jx 2 !!!,!)-^)^^ (29) 

Assuming now that the increments 5fi(i, i\) x 2 («i, i) and r(z)x 1 (z) are all elements of C% 
with /i > 1, 5fi(i, ii) x 2 (z 1; z)+r(z) x 1 (i) becomes an element of Dom(A), and inserting (12"9"|) 
into f[2"8"]) we obtain 

8z = J{m dx) = m{i) x J (i) + ii) x 2 (zi, i) — A (r(z) x 1 (z) + 5/i(z, i±) x 2 (zi, z)) , 

which is the expression (T25D of our Theorem 13.51 Thus (T25]) is a natural expression for 
J{m dx). 

Proof of Theorem Iff. 51 We will divide this proof into two steps. 

Step 1: Recalling the assumption 37 > 1, let us analyze the three terms in the right hand 
side of ( 12 5 p and show that they define an element of Q 7iQ such that 5z = ((i) x 1 (i) + f 
with 

= m(i) and f = /i(z, ii) x 2 (i l5 i) — A (r(i) x 1 (z) + 5fi(i, ii) x 2 (ii, ij) . 

Indeed, on one hand m G Cj(M. d ) and thus ( = m is of the desired form for an element of 
Q 7 , a - On the other hand, if m G Q 7 ,b, /i is assumed to be bounded and since x 2 G C2 7 (M d,d ) 
we get that n(i, i±) x 2 (z'i, i) G Cp ' ■ Along the same lines we can prove that r(i) x 1 (z) G C3 7 
and 5fi(i, i\) x 2 (z 1; i) G C3 7 . Since 37 > 1, we obtain that r(z)x 1 (z) + 5fi(i, i\) x 2 (z 1; i) G 
Dom(A) and 

A(r5x + 6fi x 2 ) G Cj 



'30 
'2 • 
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Thus we have proved that 

f = ii) x 2 («i, i) - A (r(i) x x (i) + 8fj,(i, z'i) x 2 (i i7 i)) G 

and hence that z G Q 7;0 (R). The estimate (I26p is now obtained using the same kind of 
considerations and are left to the reader for the sake of conciseness. 

Step 2: The same kind of computations as those leading to ( 12"9|) also show that 

5 (m(i) x 1 (i) + fi(i, ii) x 2 (ii, i)) = - [r(i) x 1 (i) + Sfi(i, i\) x 2 (ii, i)] . 
Hence equation (125]) can also be read as 

J{m dx) = [Id — AS] (m(i) x 1 (i) + fi(i, ii) x 2 (i l7 i)) , 



and a direct application of Corollary 12.31 yields (12 71) . which ends our proof. 



□ 



3.3. Ito-Stratonovich formula. We are now ready to obtain a change of variable for- 
mula for f(x), according to the strategy given in Section I3TT1 For this, we need to assume, 
on top of the multiplicative Hypothesis 13.11 the following geometric rule which is ([2]) in 

the case N = [±\ =2: 



-7- 

r 2 



Hypothesis 3.6. Let x be the area process defined in Hypothesis \3.1\ Then we assume 
that, for all (s,i) G Sz^t, we have 

With these assumptions in mind, our change of variable formula reads as follows: 

Proposition 3.7. Assume that x satisfies Hypothesis \3.1\ and \3.b\ Let f be a C 3 (IR d ;M) 
function. Then 

P(/(*))L = J« (V/(s) dx) = f (Vf(x u ), dx u ) Rd , (30) 

J s 

where the integral above has to be understood in the sense of Theorem \3.5l 
Proof. Consider a partition U st = {s = to < ■ ■ ■ < t n = t} of [s, t]. We have that 

n-1 



f( Xt )-f( Xs ) = J2f(xt q+l )-f(x tq ) 

q=0 

n— 1 ^ n— 1 d 

= ^9 l /(l 1 Jx t 1 A+i ( l ) + -^ 5 n, 2 /K) x M 9+ i( 2 i) x M 9+1 ( 2 2) 

q=0 9=0 ii,i2=l 

j n— 1 d 

+ d li2iJ( x tl li2i J x v 9+1 (^i) x v 9+1 fe)x t 1 9t9+i (z 3 ), (31) 



q=0 ii,i2,«3=l 

for a certain element £' liai G Invoking now Hypothesis 13.61 and Schwarz rule, 

one can express the sum \ E? ljia =i d 2 hi J(x tq ) x£ t (ii) (i 2 ) as 



d d 



\ Yl 9 ni 2 /K)[ X M 9+ l( Z l' Z 2)+X 2 ig+l (z 2 ,Z 1 )] = ^ ^ li2 /K) X V 9+1 (^1^2; 



2 

n,«2=i n,«2=i 
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So, going back to our convention on repeated indices, one can thus recast expression f[3"Tj) 
into 



n— 1 n—1 

f(x t ) - f(x s ) = J2dif(xt q )xl qtq+1 (i) + E^ 2 /W<* 9+1 (*i^) + ^ ( 32 ) 

9=0 9=0 



where _R st can be written as R st = ^ Yl^=o Pt q t q +i with 

Pt q t q+1 = df li2l JXx^) *l qtq+1 (ii) <t 9+1 (« 2 ) x^ q+1 (z 3 )- 

Furthermore, it is readily checked that for any < q < n—1 we have pt q t q+1 < C|tj+i - tj| 37 , 
with C a constant depending on / and x, owing to the fact that / is a C 3 function and x is 
continuous on [0, T\. Thus, since 37 > 1, it is easily seen that lim|n si |->o 2g=o = ®- 

Plugging this relation into fl32|) . we have proved that 



n-1 



n—1 



/(^-/(x s )=lim E^/(^ (33) 

list — ^0 



g=0 



g=0 



On the other hand, Proposition 13. 3l asserts that the decomposition of V/(x) as a weakly 
controlled path is given by 

{5Vf{x){i)) st = (SdJXx))^ = dlJ(x s )xUn) + r st (i), 

where r lies into C^ 1 Hence, using formula (12?]) . we have that 



J st {Vf{x)dx) = lim 
|n st |->o 



n—1 n—1 

E w***) <t q+1 w + E ^«/(*o x m, +1 fr» «■=») 

,g=0 g=0 



(34) 



Comparing this equality with (133|) . one obtain easily the desired Ito-Stratonovich formula. 

□ 

Remark 3.8. The following formula: 

J st (Vf(x)dx) 

E^/K)x^ 9+1 W + ^E E ^ 2 /K)<t 9+1 (H)x^ +i (z 2 ) (35) 

.9=0 9=0 11,12=1 

is also an interesting byproduct of the proof of Proposition 13.71 



lim 

|n st |-»o 



4. General Stratonovich calculus 

We will now handle the case of a weakly geometric rough path based on x G Cj(R d ) 
as defined in the introduction, and we set N = |_1 / -yj . We shall define an integration 
theory and show an Ito-Stratonovich formula for this kind of process. This being done 
along the same lines as in Section |3l we may skip some details of computations here. In 
any case, recall that we suppose the existence of a family {x 11 ; n < N} of increments in 
C2 satisfying the regularity, multiplicative and geometric properties given at Section [TJ 
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4.1. Weakly controlled processes. With respect to the case of order 2, the notion of 
controlled process is obviously obtained here by introducing more iterated integrals of the 
process x. A new kind of cascade relation is also required, which is reminiscent of the 
Heisenberg type structure of |21j . 

Definition 4.1. Let z be a process in CjiW 1 ) with l/(iV + 1) < 7 < 1/N. We say that z 
is a weakly controlled path based on x and starting at a G W 1 , if zq = a and 5z G (MJ 1 ) 
can be decomposed into 

N-l 

8z(i) = C% k, ■ ■ ■ , ik) x k fe ■■-,*!) + r°(z), (36) 

k=l 

for all (s,t) G <S2,t- In the previous formula, we assume ( k G CjiW 1 x M. dk ), and r° is a 
regular part such that r G C^iW 1 ). We also suppose that for any 1 < k < N — 2, the 
increment ( k can be further decomposed into 

N-l-k 

SC h (i,ii, ■ ■ ■ ,ik) = ^ C fc+ '(Mi> • • • *\ik+i, ■ ■ ■ ,*fe+i) + r k (i,H, . . . ,i k ), (37) 
1=1 

where the remainder r belongs to C^ - ^ 1 (W 1 x M. dk ). 

As in Section® the space of weakly controlled paths will be denoted by Q 7ia (R n ), and 
a process z G Q 7 , a (^ n ) can be considered in fact as a tuple (z, ( l , . . . , C^ -1 )- The natural 
semi-norm on <2 K ,a(K n ) is given by 

JV-l 

M[z- Q 7 ,a(K n )] = M[z\ Cj{R n )} + ^[C fe ! ^r(K n x R dk )] 

k=l 

N-l 



+M[C k ;C 1 1 (R n x R dk )] + ^AT[r fc ;Cf - fc)7 (M")], 



fc=0 



withAf\g;Cf] defined by (THJ) and Af[(; C?(V)] =sup 0<s<T |C 



s V- 



The decomposition of f(x) as a controlled process for a smooth enough function / can 
now be read as follows: 

Proposition 4.2. Let x be a path satisfying Hypothesis li.il and let f : M. d — >■ R be a 



function such that f(xo) = a, and set z = f{x). Then z G Q 7 , a ; an d it can be decomposed 
into 

N-l 

fc = ^C fc (^,---,^)x k (^,...,2i)+r , (38) 

k=l 

with 

N ~ 1 d k fir \ k 

C(h, . . . , i k ) = C,./VJ. r° = [Sf(x)] st - £ EI^^)' ( 39 ) 

k=i ' j=i 



and 



N-k-1 



1 d k+p ■ ■ fix ) p 



pi 

p=l q=l 
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where we recall our convention ( fFTj) for Q^...^ /. Furthermore, 

N[z\ Q 7 , a } < c f>T (1 + A^fx; C?(R d )]) . (41) 

Proof. The algebraic part of the assertion is obtained by combining a simple Taylor ex- 
pansion and our geometric assumption (j2J). Indeed, Taylor's expansion directly yields 

N — l rjfc /• / \ k 

[Sf(x)] st = D st + r° st , with D st = J2 ki ( 42 ) 

fe=l " j=l 

where r° G C^ 7 . Moreover, appealing to (J2J), we have 

JV-l 1 d k 

= El E #,V0O (43) 

fc=l ii,—>*fc=l J=l 
JV-l . d 

= E^i E ti/WE^)-''^)). ( 44 ) 
fe=l ' h,...,i k =i o-es fc 

and invoking the symmetry properties for the derivatives of / we obtain 

JV-l . c2 

Drf = E&i E E a t w -vi,/( :I ') xk (Mi).-.Mfc)) 
fc=i ii,...,j fe =i o-eSfc 

JV-l d 

= E E C ,/'(•'•>) xk «'i (45) 

fc=l U,...,j fe = l 

Going back to our convention on repeated indices, we end up with 

JV-l 



which is the announced formula (l38|) . Further expansions of the coefficients 9^ il f(x), 
leading to a relation of type (}3"Tj) . are performed in the same way and are left to the reader 
for sake of conciseness. 

In order to give an estimate for Af[z; Q 7ja (IR n )], one has of course to establish bounds for 
M[z\ CjiW 1 )], A/"[C fc ; C7(E d )], Af[C k ;C^°(R d )] and A/> fc ;C 2 27 ]. These estimates are readily 
obtained from the expressions in decomposition ( 15§|) . and details are left to the reader. 
The analytic bound (|4ip is also obtained in a straightforward manner. 

□ 

Remark 4.3. As in the case n = 2, we point out that the algebraic conclusion of this 
proposition is still true for a / G C N (M. d ), without boundedness restrictions. However, 
inequality (14T1) would take a different form, since the multiplicative constants depend on 
the derivatives of / composed with x. 

4.2. Integration of weakly controlled paths. The formula which defines the integral 
of a controlled process with respect to x is now defined similarly to the one in Theorem [321 
in spite of the roughness of x. 
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Theorem 4.4. For a given 7 > with [1/tJ = N (that is, 1/(N + 1) < 7 < l/N), let x 
be a process satisfying Hypothesis Furthermore, let m G Q 7i f,(]R d ) with decomposition 
nio = b G M. d and 

N-1 

8m{i) = h,---, ik) x k (4, ...,h) + r°(i), (46) 

k=l 

where the increments n k satisfy the further assumptions of Definition ^..1 Define z by 
zq = a G R and 

N-1 

5z = m{i) x 1 ^) + *i> • • • > «'fc)x k+1 (4, . . . , ii, i) (47) 

k=l 

(N-2 
^2 rh (h H, ■ ■ ■ , 4)x k+1 (i fc , . . . , i u i) + 5^ N ~ l (i, h,..., iN-Jy^iiN-i, ■ ■ ■ , h, i) 
k=0 



Finally, set 



J st (mdx)= / (m u , dx u ) Rd = (5z) st . 



Then: 

(1) z is well-defined as an element of Q Kj0 (R), and coincides with the Riemann integral 
of z with respect to x whenever these two functions are smooth. 

(2) The semi-norm of z in Q K a (M) can be estimated as 

M[z- Q 7 ,„(E)] < c x (1 + Af[m; Q.y, b (R d )}) , (48) 

for a positive constant c x which can be bounded as c x < c^^ x A^[x k ; C 1 ^ 1 ], where c stands 
for a universal constant. 

(3) It holds 

n-1 N-1 

J st {mdx) = lim V \ m tq (i) ^ tq+1 (i) + ^ ^(i, k,..., i k ) xj^fo, • • • ,M) (49) 
|nstK %=o L k=i J 

for any < s < t <T, where the limit is taken over all partitions U s t = {s = to, . . . , t n = 
t} of[s,t], as the mesh of the partition goes to zero. 

Proof. Relying on what has been done at Section [31 we mainly derive here the expres- 
sion f l4"Tj) for Jimdx). Once this expression is obtained, the other estimates follow like 
in Theorem 13.5} except for the higher number of terms which have to be taken care of. 

Hence let us assume for the moment that m and x are smooth functions, and try to 
define J{m dx) in an appropriate way for generalizations to rougher cases: one can write, 
using decomposition f l4"6]) . 



J{mdx) = m(i) x 1 (i) + J (Sm(i) dx) (50) 

N-1 

m{i)^({) + ^/i fc (z,2 1 ,...,z fc )x k+1 (2 fc ,...,z 1 ,z) + J(r°dx). (51) 

k=l 
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Hence, like for equation fl2UJ), one can deduce that 



N-l 

.k/- ■ ■ \„k+l, 



5(J(r°dx)) = -5m(i)x 1 (i) - ^ 5[i k [i, i u . . . , i k )yi" 1 -[i k , . . . , zi, i) 

k=l 

N-l 

+ ^ A**(*\ *i> ■ ■ ■ , 4)5x k+1 (i fc , . . . , i t , i). 



k=l 



We now plug relation ()46|) for 5m, relation f l37|) for <5/i fc and the multiplicative relation §1§ 
for 5x k+1 into the latter equation. This yields 



N-l 



S(J( r ° dx)) = - h, ■ ■ ■ , ik) x k (^, . . . , h) x^i) - r°(i) x 1 ^) - M 

fc=i 

AT-2 

- cfyt N-1 (i, ijv-i) x N (zV_i, ■ ■ ■ - ^ r fc (z, ii, . . . , x k+1 (i fc , . . . , z'i, i) 



fc=i 

AT-1 fe 

+ fi k (i, i l7 ...,i k )^2 x '(^> • • • > xk+1_1 (*fc-i) — , ii, i), (52) 

fe=i ;=i 

where 

JV-2 N-l-k 

M =Y{ ^2 /-^(Mi, • • • ' ifc +') x '(^+«, • • • ,ik+i)j x k+1 (z fc , . . . ,h,i). 

k=l 1=1 

Moreover, a simple change of index allows to write 

N-l q-l 

M = 0*(*> • • • , i q ) x '(v • • • ' Vi+i) x q+1-1 (Vi> 

q=2 1=1 

and hence (152]) simplifies into 

N-2 

= ~ *!'•••» ^)x k+1 (i fc , . . . , h, i) - 8n N ~ l {i, ii, . . . , i A r_i)x N (i 7V _i, — , i x , i). 

fc=0 

(53) 

It is now readily checked that the operator A can be applied to the latter increment when- 
ever x G CI and generates a weakly geometric rough path. Putting together relations (|50j) 
and (l53j) we thus end up with expression (|47|) for the integral J(mdx). 

The analytic bounds are now a matter of standard calculations, and are left to the 
reader for sake of conciseness. 

□ 



4.3. Ito-Stratonovich formula. Now that we know how to define integrals of controlled 
processes with respect to x, our change of variable formula for f(x) is obtained quite in 
the same way as in the second order setting. The formula can then be read as follows: 
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Theorem 4.5. For a given 7 > with |1/tJ = N, let x be a process satisfying the 
regularity, multiplicative and geometric hypotheses of Section [3 Let f be a C N+1 (R d ; M) 
function. Then 

W(x))} st = Jst(Vf(x)dx)=[ (Vf(x u ),dx u ) Rd , (54) 



where the integral above has to be understood in the sense of Theorem 4-4 Moreover, 

71-1 N-1 

J st (V/(x) dx) = Jim J] \dj(x tq ) xl tg+1 (i) + 9t +1 ,JM ■ ■ ■ '«'i»0 

list — HJ L 

<3=0 fc=l 

n-1 N-1 



n-1 JV-1 . 

n U t f, S [ 9 </( z O x M 9+ i( 2 ) + S ^[Cii^K) x V 9+ i(^) • ■ ■ x M, + i(^) X V 9+1 < 



" q=0 " fc=l 



(55) 



/or any < s < t < T, where the limit is taken over all partitions U s t = {s = to, . . . , t n = 
t} of [s,t], as the mesh of the partition goes to zero. 

Proof. The proof goes exactly along the same lines as for Proposition 13.71 The first 
expression for J(¥ f(x)dx) st given in (|55j) is a consequence of the decomposition as a 
weakly controlled path of V/(x) and Theorem 14.41 The second one follows from the first 
one by using Schwarz rule and the geometric property (121) . 

On the other hand, applying Taylor's formula up to order iV to the decomposition 



n-1 

f{x t ) - f{x s ) = f{xt q+1 ) - f(x tqJ 

q=Q 



comparing it with (135]) and using that (N + 1)7 > 1, one obtain easily the Stratonovich 
type formula. 

□ 



5. SKOROHOD TYPE FORMULA VIA MALLIAVIN CALCULUS 

We take now a completely different direction in our considerations: the pointwise point 
of view which had been adopted previously is abandoned in this section, and we try 
to construct an integral with respect to a (Gaussian) process x by means of stochastic 
analysis tools. We then prove that for any < s < t < 00, the function l[ s>t )V/(x) is in 
the domain of an extended divergence operator with respect to x, and prove an associated 
Skorohod type formula. As we shall see, this mainly stems from an extension of j2] to the 
d-dimensional case, which is allowed thanks to the symmetries of Vf(x). 

5.1. Preliminaries on Gaussian processes. >From now on, we specialize our setting 
to a centered Gaussian process x = (x(l), . . . , x(d)) with i.i.d coordinates, and covariance 
function 

R st :=E[x 8 (l)z t (l)], and R t := E[\x t (l)\ 2 } = R tt , s,te[0,T\. (56) 
We will add later some hypotheses on these functions. We can also assume that Xo(j) = 0. 
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The Gaussian integration theory is based on a completion (in L 2 (Q)) of elementary 
integrals with respect to x, which can be summarized as follows (see [25] for more details): 
consider the space of d- dimensional elementary functions 

nj — l 

s = {/ = (h, ■ ■ ■ , fd)\ fj = 4\i,i +l) , o = t < 4 < ■ ■ ■ < < ti, = t, 

for j = 1, . . . ,d\ . 

For any element / in S, we define the integral of first order of / with respect to x as 

j=l t=0 

For 9 : R ->• K, and j G {1, . . . , d}, denote by 00 the funct ion with values in M. d having; 
all the coordinates equal to zero, except the j th coordinate which is equal to 9. It is readily 
seen that 

E [ / i( 1 [o, a )) / i( 1 |o!t))] = 1 (j=k)Rst- 
So, we can define for some indicator functions of S the following symmetric and semi- 
definite form 

h b'l i [ fc l \ _ i p 
\- L [o,s)' ^[ofyiS — - L 0=fc)-n-st, 

and extend it to all elements of S by linearity. If we identify two functions / and g in S 
when (/ — g, f — g)s = 0, then (•, becomes an inner product on S (actually, on the 
quotient space obtained by this identification). Therefore, for / and g in S we have that 

E[I 1 (f)h(g)} = (f,g)s 

and Ji defines an isometric map from S, endowed with the inner product (•, •)$ into a 
subspace of L 2 (Q). This map can be extended in the standard way to an isometric map, 
denoted also as Ii, from a real Hilbert space that we will denote by % into a closed 
subspace of L 2 (Q). >From now on, denote the inner product of this extended isometry 
by {■■,■)%■ We will assume that H is a separable Hilbert space (which is satisfied whenever 
R s t is continuous). 

Let {ei, C2-, ■ ■ ■ } be an orthonormal basis of % and let ® denote the symmetric tensor 
product. Then 

U = /u,-,^e n ® ■ ■ - ®e in , f ilt ... )in e M (57) 

finite 

is an element of '}i® n with the Hilbert norm 

||/n||^/<g)n — / ] \fii,— ,in\ 2 ■ (58) 
finite 

Moreover, %® n is the completion of all the elements like (1571) with respect to the norm (l58]h 

For an element f n G %® n , the multiple Ito integral of order n is well-defined. First, any 
element of the form given by fl57|) can be rewritten as 

fn = J2fn---iA kl ®---®€ m > ( 59 ) 

finite 
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where the ji, . . . , j m are different and k\ + ■ ■ ■ + k m — n. Then, if f n G "H®" is given under 
the form fl59|) . define its multiple integral as: 



Wn) = J2fn,-^H kl (I 1 (e n ))..-H k JI 1 (e 1 J), (60) 

finite 

where H k denotes the k-th normalized Hermite polynomial given by 

Hk(x) = = V . iz^l - x *-*. 

v ; v 1 dx k ^ 23j\(k-2j)\ 

It holds that the multiple integrals of different order are orthogonal and that 



E\I n (fn)\ 2 = n\ \\f n \\ 2 



Y II 



This last isometric property allows to extend the multiple integral for a general f n G FL 
by L 2 (Q) convergence (notice once again that this kind of closure is different in spirit 
from the pathwise convergences considered at Sections |3] and H]). Finally, one can define 
the integral of f n G 1-L® n by putting I n (f n ) := I n (fn), where /„ G V® n denotes the 
symmetrized version of f n . Moreover, the chaos expansion theorem states that any square 
integrable random variable F G L 2 (Q,Q, P), where Q is the a- field generated by x, can 
be written as 

oo oo 

F = Y, I Mn) with E[F 2 ] = 5>!||/„||^. (61) 

n=0 n=0 

We will introduce now the (iterated) derivative and divergence operators of the Malli- 
avin calculus. We denote by Cf°(lR n ) the set of infinitely continuously different iable func- 
tions / : MJ 1 — > K such that / and all its partial derivatives have polynomial growth. Let 
S denote the class of smooth random variables of the form 

F = f{I 1 (h 1 ),...,I 1 (h n )), (62) 

where / G C^°(K n ), h\, . . . ,h n are in "H, and n > 1. The derivative of a smooth random 
variable F G S of the form ( 162]) is the H- valued random variable given by 

n 

DF = J2 %f (11(h), h(K))K (63) 
i=i 

where di denotes as usual One can also define for h G T-L and F G S the derivative of 
F in the direction of h as D^F = (F, h)u. 

The iteration of the operator D is defined in such a way that for a smooth random 
variable F G S the iterated derivative D k F is a random variable with values in T-L® k . We 
also consider for h k G r H® k the fc-th derivative of F in the direction of h k defined as 



D k hk F = {D k F,h k ) 



k. 



Let us fix now a notation for the domain of the iterated derivative D . for every p > 1 
and any natural number k > 1 we introduce the seminorm on S given by 



E(\F\*>)+J2H\\D j F\\ p mj ) 
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It is well-known that the operator D k is closable from S into L P (Q; T-L® k ). We will denote 
by 3 k,p the completion of the family of smooth random variables S with respect to the 
norm || • \\k, P ■ We will also refer the space D fc ' 2 as the domain of the operator D k and 
denote it by DomD fc . If F has the chaotic representation (1611) . we have that 



E(||£> fc F|& w ) = - 1) • • • (n - k+ 1) n\\\f n f u&n 

n=k 

and a useful characterization of DomD fc is the following: F G DomD fc if and only if 

oo 



n=l 

We will denote by 5° the adjoint of the operator D (this operator is also referred as 
the divergence operator) and more generally, we denote by 5 ofc the adjoint of D k . The 
operator 5 ofc is closed and its domain, denoted by Dom 6° k , is the set of "H® fc - valued square 
integrable random variables u G L 2 {VL\'H® k ) such that 

\E((D k F,u) H m\ <C||F|| 2 , 

for all F G Dom D k , where C is some constant depending on u. Moreover, for u G Dom S 0k , 
S 0k (u) is the element of L 2 (Q) characterized by the duality relationship: 

E(F5 ok (u)) = E((D k F, u) mk ), (64) 

for any F G DomZ) fc . For u G Dom 5°, the random variable 5°(u) is usually called 
Skorohod integral of u, because it coincides with the usual integral of u with respect to x 
for a large class of elementary processes u (see [25] for further details). 

5.2. An operator associated to x. Along this section we will consider a <i-dimensional 
continuous process satisfying the following set of assumptions: 

Hypothesis 5.1. The process x = (x(l), . . . , x(d)) is a centered Gaussian process with 
i.i.d. coordinates. Letting R st and R t being defined as in l[5h]) . we suppose that those two 
functions are continuous and the following two conditions hold: 

(1) The variance function R t := R tt is differentiable at any point t G (0,T) and 
satisfies that 

\R'.\ dt < oo. 



(2) The first partial derivative d s R s t of R st is well-defined a.e. on [0,T] 2 and verifies 

/ / \d s R sy \ dsdy < oo. (65) 
Jo Jo 

We will try now to identify a useful operator for our future Gaussian computations. 

Let <p = (y?(l), . . • , <p(d)) G (T*T) d , where T>^ is the space of C°° functions with compact 
support contained in (0,T). We have that (see for instance [18], where the 1-dimensional 
case is considered) that ip G % and that 

J i(y?) = - / fa, v's) ds i 
Jo 
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where (-, ■) denotes the ordinary Euclidean product in M. d . Moreover, (T>T) d is a dense 
subset of %. From now on, we use also the notation x(f) for ii(/). 

Given a function h : [0,T] — > JR., recall that denotes the function with values in M. d 
in which all the coordinates except the j-th one are equal to and the j-th coordinate 
equals to h. Therefore, for G XV and < a < b < T, we have that 

(l [atb f,^% = E[J 1 (l [a>6) M)J 1 G9W)]=-Ef( i F 6 (Z)- a ; a (0) [ T (x t ,[(3%)dt 

Jo 

= - 10=1) J q (Rbt - Rat) P' t dt = - l (j=l) jf Qf d s R st ds^j P' t dt 

= -l U= i) J n d s R sy (3' y dy^ds. (66) 

We will consider the first iterated integral appearing on the right hand side of f lB"6"|) as 
a linear operator defined on XV- That is, we consider for s G [0, T] and G XV, the 
following function: 

A0(s) =: - / d s R sy (3' y dy. 
Jo 

We will suppose from now on that the following hypothesis holds. 
Hypothesis 5.2. For any (3 G XV, A/3 G X 2 ([0,T]). 

Remark 5.3. Condition (I6"5j) on i2 st , stated in Hypothesis l5.lt implies that A/3 belongs to 
L 1 ([0,T]) whenever /3 G XV- We have imposed the additional condition Af3 G L 2 ([0,T]) 
in order to guarantee the integrability of many terms appearing in the sequel. Although 
one can weaken Hypothesis 15. 2\ this would complicate some of the next statements. We 
have thus chosen to impose it for the sake of simplicity. 

Example 5.4. Hypothesis 15.21 is satisfied by the fractional Brownian motion. In fact, 
A/3(s) = - J T ^R sy P\y)dy = - j\(s 2H - l -\s-y\ 2H ~h\gn{s-y)) 0'(y)dy 

= [ H(\s-y\ 2H - 1 sign(s-y))/3'(y)dy, 
Jo 

because /3 G XV- And from this, it is easily seen that A/3 G L°°([0,T]) if (3 G XV- 

Let us now relate our operator A to the inner product in %: equation (166|) tells us that 
for any elementary function g = (g(l), . . . , g{d)) G S and G XV we have that 

(0, g(l)M) H = l u=l) f g s (l)A0(s)ds. 



Since for ip G (XV) d , g G S, we have </? = Y^j=ifU)^ ano - S 1 = Sf=i°(0^5 we obtain 
that 

d „t pT 

( < P>9)n = yZ 9s(j)ALp(j)(s)ds= (g s , A<p(s)) ds, (67) 
J=1 </o Jo 

where we use the notation Aip = (A<^(1), . . . , Aip(d)). Extending this last relation by 
continuity, the following useful representation for the inner product in "H is readily ob- 
tained: 
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Lemma 5.5. For any g G Ji H (L 2 ([0,T])) d and tp G T>t> one can write 

(<P,9)h = (9s, A.<p(s)) ds, (68) 
Jo 

where (•, •) stands for the inner product in M. d . 



Going back to our example 15.41 notice that expression (168J) is similar to the following 
one pointed out in |2] for the one- dimensional fractional Brownian motion with Hurst 
parameter H < 1/2: 



(<P,9) H = <& [ g(s)r> a + -D a Ms)ds, 
Jo 



where a = | — H; D" and D" are the Marchaud fractional derivatives (see |28j for more 
details about these objects), and c# is a certain positive constant. 

5.3. Extended divergence operator. Let us take up here the notations of Section 
15.11 Having noticed that ffim gives rise to an operator D"D" which is a particular case 
of our operator A (see Example I5.10p . one can naturally try to define an extension of 
the operator 5° using similar arguments to those of [2J. The idea is to consider first 
u G Dom5° n (L 2 (f2 x [0,T])) and F = H n (x{ipj) where H n is the n-th normalized 
Hermite polynomial, and tp G (V T ) d . Since (L 2 (fi x [0,T])) d = L 2 (Q; L 2 ([0, T]; R d )) 
and Dom<5° C L 2 (fi;"H), we have that u G (L 2 (fi x [0,T])) d n"H almost surely. Moreover, 
DH n -i(x(<p)) = H n _i(x{ip))(p G (T> T ) d , a.s.. So, using ( 168]) . the usual duality relationship 
between D and 5° can be written in the following way: 

E [S*(u)H n (x{<p))] = E [{«, DH n (x(<p))) n ] = E [H n ^(x(if)) (u,<p) n ] 



E 



H n -\ (#(<£>) ) / (u a , A<p(s)) ds] = [ (E [/f n _ 1 (x(^))«J , A<f(s)) ds, (69) 
Jo J Jo 

and this motivates the following definition. 

Definition 5.6. We say that u G Dom*5° if u G (L 2 (f2 x [0,T])) d and there exists an 
element of L 2 (Q), that will be denoted by 5°(u), such that for any (p G (XV) an d an V 
n > 0, the following is satisfied: 

E [5\u)H n (x{v))] = [ T (E [H n ^(x^))u s ] , Acp(s)) ds. (70) 

Jo 

Remark 5.7. Since the linear span of the set {H n ((p) : n > 0, tp G {VtY} is dense in 
L 2 (Q), the element 5°(u), if it exists, is uniquely defined. 

Remark 5.8. One can easily see from our definition of the extended divergence that it is 
a closed operator in the following sense: if {w fc }fceN C Dom*5° and satisfies (1) u k — > u in 
(L 2 (tt x [0,T})) d and (2) 5°{u k ) -> X in L 2 (fi), then u G Dom*<5° and 5°(u) = X. 

We show in the following proposition that the extended operator 5° defined above is 
actually an extension of the divergence operator of the Malliavin calculus. 

Proposition 5.9. The domain Dom*S° is an extension of Dom5° in the following sense: 

Domb° n (L 2 (n x [0, T])) d = Dom*5 <> n L 2 (fi; U). 
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Furthermore, the extended operator 5° restricted to Dom5° fl (L 2 (f2 x [0,T])) coincides 
with the standard divergence operator. 

Proof. Uu G Dom^n(L 2 (fi x [0, T])) d then u G (L 2 ([0, T])fnl-L almost surely. Thus ([68]) 
can be applied to u and ( 169]) holds true for (the standard divergence operator). This 
proves that Dom^ n {L 2 (Q x [0, T])) d C Dom*^ n L 2 (fi; H) and that 5° is an extension 
of the standard divergence operator on Domtf R (L 2 (VL x [0,T])) . 

To see the other inclusion, take u G Dom*5° fl L 2 {Vl\FL). By our Definition 15.61 of 
Dom* 5°, u belongs also to (L 2 (f2 x [0,T])) . We will show that u G Dom<5°. First, we 
will prove that the element 5°(u) defined by the equality (170]) satisfies, for any ip G {T>T) d 
and any n > 0, that 

E [8*(u)H n (x((p))] =E[{u,DH n (x(<p))) H ]. (71) 

Indeed, since u G (L 2 ([0, T])) d D"H a.s. by assumption, we can apply again identity (|68|) 
and so, 

(m , (x(v?)) ) H = H n -! (x(ip)) / (u s , Ay?(s)) ds. 

■/ o 

Hence, using Fubini's theorem and (I70p we end up with 

E(u,DH n (x(ip))) = f T (E \H n -i{x{ip)) u s ] ,A v (s))ds = E[S%u)H n (x(ip))], 

Jo 

which is exactly (1711) . 

By using density arguments (the linear space generated by the elements of the form 
H n (x(ip)^, with if G (T>T) d , n > 0, is dense in DomD) we obtain that 

E[(u,DF) H ]=E[5%u) F] 

for any F G DomD, and this finishes the proof. 

□ 

Example 5.10. Go back to our ffim Example 15.41 and let us compare the extended diver- 
gence operator introduced above with the one defined in [2]. First of all, we must point 
out that in [2], the (standard) divergence operator is presented in a more general setting 
than ours: the divergence can belong to any L P (Q), for p > 1. In our paper, we will only 
consider this divergence over L 2 spaces for sake of conciseness. 

According to the computations carried out in [16J (see identity (5.30) of that work), for 
any element ip in the space of test functions T>t one has 

c|D^D^( S ) = / H\s - y\ 2H -hign{s - y)^(y) dy. 
Jo 

On the other hand, we have already seen at Example 15.41 that Atp(s) = H \s — 
y\ 2H ~ 1 sign(s — y)ifj'(y)dy. That is, on T>t, we have the following identity of opera- 
tors: c^D"D" = A. Moreover, these operators can be extended (and coincide) by 
density arguments to I°:(Sh) (see [2] for the definition of this space). Finally, in this 
case, FL = /°(L 2 ([0,T]) is a subset of L 2 ([0,T]). Using these observations, it is readily 
checked that the extended divergence operator defined above coincides with the extended 
divergence given in [2J, restricted to L 2 spaces. 
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5.4. Change of variable formula for Skorohod integrals. We can now turn to the 
main aim of this section, namely the proof of a change of variable formula for f(x) based 
on our extended divergence operator 5°. Interestingly enough, this will be achieved under 
some non restrictive exponential growth conditions on /. 

Definition 5.11. We will say that a function f : M. d — > R satisfies the growth condition 
(GC) if there exist positive constants C and X such that 

X< — -, and \f(x)\< Ce x lxl " for all x eM d . (72) 

Ad max te[0 ,T] Rt 

Notice that max te [o 5 T] Rt = niax tg [o,T] -^[l^tl 2 ]- Thus the growth condition above implies 
that 



E[max \ f(x t )\ r ] < C r E^e rXmaXt ^ lxtl 



te[o,T] 

and this last expectation is finite (see, for instance |23| . Corollary 5.4.6) if and only if 

1 1 

rX < 



2max te [ ,T] E(\x r \ 2 ) 2d max tG [ ,r] Rt 
So, if condition (1721 is satisfied, there exists r > 2 such that 



E 



max \f(xt)\ r < oo. (73) 
te[o,T] 



With these preliminaries in hand, we first state a Skorohod type change of variable 
formula for a very regular function /. 

Proposition 5.12. Let f £ C°°(IR d ) such that f and all its derivatives satisfy the growth 
condition (GC) (with possibly different X's and C's). Then, for any < s < t < T, 

l M (-)V/(x.)£ J Dom*^ 

and 

5" [1 M (-) V/(s.)] = f(x t ) - f{x.) -\j Af(x p ) R p dp. 

Proof. Since / and all its derivatives satisfy growth condition (GC) , the process l[ s ,t) V/(x) 
is an element of (L 2 (fi x [0,T])) d and we also have 



f(x t ) - f(x s ) -\J Af(x p ) R' p dp e L 2 (Q). 



So, we only need to show that for any n > and any ip £ (T>t) the following equality is 
satisfied: 



E 



f(x t ) - f(x s ) -\\ Af(x p ) R p dp ) H n (x(<p)) 



[E [H^x^)) V/(x p )] , A<p(p)) dp. (74) 



The proof of this fact is similar to that of [21 Lemma 4.3], although some technical 
complications arise from the fact that here we deal with the multidimensional case. 
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Consider thus the Gaussian kernel 



_ a 

p(a, y) = (2-na) 2 exp 



i\yl 

2 a 



for a > 0, y G 



(75) 



It is a well-known fact that d a p = \Ap. Moreover, E = J Rd p(Rt,y) g(y) dy for any 

regular function g : M d — > M. such that g and all its derivatives satisfy (GC) . Using these 
identities, we can perform the following computations: 



d_ 

dt 



(i 

dt 



p(Rt,y)g(y) dy 



d_ 
~d~o 



p(Rt,y) R' t g(y) dy 



l -R' t ( Ap(R t ,y)g(y)dy=]-R' t [ p(R t ,y) AG(y)dy 



-R' t E[Ag(x t )]. 



(76) 



This shows that the function ^E[g(xt)] is defined in all t G (0, T) and is integrable on 
[0,T] . As a consequence, Efg^)] is absolutely continuous. Using this fact and identity 
(176)) . we can now prove (l74j) when n — 0. Indeed, observe that in this case H (x) = 1 
and, by definition, H~i(x) = 0. Hence, the right-hand side of (174|) is equal to while the 
left-hand side gives: 



E 



f( Xt )-f(x s )-- / Af(x p )R'dp 



E[A/(x„)] Bf.dp, 



and this last quantity vanishes due to (I76p . 

Let now n > 1. Define for j G {1, . . . , d} and t G [0, T], 

rp 



GUp) 



A<p(j)(s) ds 



;iM W .vCf) M >„- 



(77) 



Clearly, G^, is absolutely continuous and (G^,) = A.<p(j) (a.e). Moreover, for a regu- 
lar function g satisfying (GC) together with all its derivatives and for any multiindex 
(ji, ■ ■ • , 3n) e {1, . . . , d} n we have 



j- p (E [g{x p )\ G»(p) ■ ~G*{p)) = i E [A^)] # p G*(p) • ■ ■ G*(p) 

+ E[g(x p )}J2^<fiUr)(p) 

r=l 



(7* 



where we have used (176]) . Recall now our convention (121 p . allowing to write <9™ - n / for 



/, and set M 31 = E[^.„ in /(x p )] FGLi^W for p G By integrating (J78 



dyj 1 --dy jn ■ 

from s to £ and taking o = <9™ ■ f we obtain that 



M 3l...J« 



Ml 1 " 3 " = - 
2 



- jT (e [dl.. jn f(x p )] (f><p(; r )(p) II ^(P))j dp. 



I: lytr 
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Summing these expressions over all the multiindices . . . , j n ) G {1, . . . , d} n and owing 
to the fact that 

/t n 
EK,,„/W](EM>)n G 5W)^ 

n ,-,Jn - r=l hl^r 

rt n-l 

/ E[d£^f(x p )]Y[G*(p)Atp(j n )(p)dp, 



= n 

31,-,3« ' 

we end up with an expression of the form 



1=1 



1 f % n 

e w 1 "' jn - ■"'"] = 2 e / E [ Ad i..inf(xp)] KH G i(p) d p 

1,-Jn 31,—>3n S '=1 

E II^^) M>Un)(p)dp. (79) 



31,— ,3n 



It should be observed at this point that, as in identity flo3]) . the symmetries of the partial 
derivatives of / play a crucial role in the proof of the current proposition. This symmetry 
property appears precisely in the computations above. 

We will see now how to obtain the desired identity (1741) from (1?T)|) . Indeed, it is a 
well known fact (see |25| again) that H n _ 1 (x((p)*) (p G Dom5° and <5° \H n _i(x{<pj) tp\ = 
nH n (x(ip)^}. Using these last two facts, the duality relationship between D and 5° and 
the definition (I77p of we have that, for g satisfying (GC) as well as its derivatives, 



E [H n (x(tp))g(x t 



-E(H n ^(x(^))^, Dg(x t )) H 
1 



Ie 

n 



E (x(<p))<p, l [0tt )Vg(x t )) n 

d rt 

H n -i(x(ip)) E / djg(x t ) Aip{j){p)dp 



3=1 



1 ^ 



3=1 



Iterating this procedure n times, one ends up with the identity 

E [H n (x(v))g(x t )] = I E E ■ ■ 



3l,-3Ti 



As an application of this general calculation, we can deduce the following equalities: 

E[tf n (x(^))/(* t )] = iE^^WW-^O ( 8 °) 



31,-Jn 



E [^(xfoOj/C*.)] = ^E[^ S )]GJ( S )...GJ( S ) (8i; 



31,-Jn 



E [H n (x(<p))Af(x p )} = ^Eft^MW-W 



f82] 



31,-3* 
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and 

E [H^x^djfM] = -^-y Yl E l^h^-dfM] G £(p) ■ • -G^ip). (83) 
Substituting now (I80|) - (l83|) in (179|) . we obtain 
n\E [H n (x(tp))f(x t )] -n!E [tf n (s(^))/(x a )] = ^n! J E [i/ n (i(^)) A/(s„)] i^tZp 



+ n (n - 1)! f E [^i(^)) 0in/(*p)] dp, 



jn = l 

and this is actually equality (j74|) . The proof is now finished. 



□ 



Since the Skorohod divergence operator is closable, we can now generalize our change 
of variable formula: 

Theorem 5.13. The conclusions of Proposition \5.1S\ still hold true whenever f is an 
element of C 2 (M. d ) such that f and its partial derivatives up to second order verify the 
growth condition (GC) . 

Proof. Let A be the constant appearing in the growth condition (GC) . Given k > 2A, 
denote by Pk(y) = piziV) the Gaussian kernel defined in (1751) and introduce fk(y) = 
(/ *Pk)(y) where, as usual, * denotes the convolution product. 
We first claim that there exist kn G N, C > and A' satisfying A < A' < -ri 5-, 

u ' 4d max tg [ ,T] Rt ' 

such that 



sup 

k>kg 

Indeed, condition (GC) easily yields 

d/2 



\fk{y)\<C'e x '\y\\ (84) 



\ a/z « 2 

f / \}{y-z)\e- kJ ^dz 

< (If c L*"'** - 4 (^// (K -* ,a A < 



On the other hand, 




2tt 

and lim^on r^rrr = A. Hence, given A' G (A, -r-; 5-), there exists L 6 N such that 

for any k > ko the following inequalities are satisfied: 



A k 1 
A < — < A < 



k — 2A 4 c? max tG [ 0i r] Rt 

Our claim (184)) is now easily deduced. 



STOCHASTIC CALCULUS FOR GAUSSIAN PROCESSES 



29 



Notice that (I8~4"P means that for k > fco, fk also satisfies the growth condition (GC) (with 
C and A' substituting C and A, respectively). Moreover, we have that 



E 



sup sup \fk(x p )\ 2 

P6[0,T] k>k 



< oo. 



Thanks to this inequality, as well as similar ones involving the derivatives of /, one can 
easily see that: 

(1) fk{x a ) -> f(x a ) and f k (x t ) -> f(x t ) in L 2 (Q), 

(2) £ Af k (x p ) R' p dp £ A/(x„) R' p dp in L 2 (fi) and 

(3) l [s , t) Vfk(x.) -+ l [s ,t)V/(x.) in (L 2 (fi x [0,T])) d . 

The result is finally obtained by applying Proposition 15.121 and the closeness of the ex- 
tended operator 5° alluded to at Remark 15.81 

□ 



6. Representation of the Skorohod integral 

Up to now, we have given two unrelated change of variable formulas for f(x): one 
based on pathwise considerations (Theorem 14. 5p and the other one by means of Malliavin 
calculus (Theorem I5.13p . We propose now to make a link between the two formulas and 
integrals by means of Riemann sums. 

Namely, let x be a process generating a rough path of order N. We have seen at 
equation (155]) that the Stratonovich integral J~ s t(V f(x)dx) is given by limin st i->.o S nst , 
where 

n-l N 

sYlat ■= E E ^Cu/K) x U +1 (*i) x U +1 fe) • ■■xi tg+1 (i k ). 

q=0 fc=l 

In a Gaussian setting, it is thus natural to think that a natural candidate for the Skorohod 
integral 5°(V/(x)) is also given by lim|n st |->.o S UBt '°, with 

n N 1 

S Uat '° ^EE^t^/W^^JiOo- -o^UJik), (85) 

g=0 fc=l 

where o denotes the Wick product. We shall see that this is indeed the case, with the 
following strategy: 

(i) One should thus first check that linin st S Uet,<> exists. In order to check this convergence, 
we shall use extensively Wick calculus, in order to write 

dt..ij{xt q ) «< t , +1 (ii) o • • ■ ox^ +i (z fc ) = d^ J\x tq ) x^ 9+1 (zi) • • -x^ +i (z fc ) + p tgjtg+1 , 

where p is a certain correction increment which can be computed explicitly. Plugging this 
relation into (I85p . we obtain 

n-l 

s n 3t ,o = s n 3t + Y^p tq , q+1 . (86) 

9=0 

(ii) Manipulating the exact expression of the remainder p, we will be able to prove that 
lini|n st |^o YlqZi Pt q ,t q+1 = ~\ Is Af(x v ) R' v dv. Hence, going back to QBE) and invoking the 
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fact that S Ust converges to J' s t{Vf(x) dx), we obtain 

lim S u ^ = J st (Vf(x) dx) -\f Af(x v ) R' v dv = [8f(x)] st - \ f Af(x v ) R' v dv. 
|n si |->-o 2 J s 2 J s 

This gives both the convergence of S nsU<> and an Ito-Skorohod formula of the form: 

[5f(x)} st = lim S u ^ + I f Af(x v ) R' v dv. 
|n st |-s-o 2 J s 

(iii) Putting together this last equality and Theorem 15.131 it can be deduced that under 
Hypotheses 15.11 and 15 . 2 \ the limit of S st ' coincides with the Skorohod integral 5°(Vf(x)). 
This gives our link relating the Stratonovich integral J{W f{x) dx) and the Skorohod 
integral 5°(V/(x)). 

This relatively straightforward strategy being set, we turn now to the technical details of 
its realization. To this end, the main issue is obviously the computation of the corrections 
between Wick and ordinary products in sums like S net,<> . We thus start by recalling some 
basic facts of Wick computations. 

6.1. Notions of Wick calculus. We present here the notions of Wick calculus needed 
later on, basically following |17| . We also use extensively the notations introduced in 
Section 15.11 

One way to introduce Wick products on a Wiener space is to impose the relation 

In(fn) O Im{9m) = In+rn[fn®9rn) 

for any /„ G 1i® n and g m G "H®' m , where the multiple integrals I n (f n ) and I m (g m ) are 
defined by . If F = I n (f n ) and G = £^ 2 =1 I m (g m ), we define FoGby 

Ni N 2 

F OG = ^^2 I n+m(fn®9m)- 
n=l m=l 

By a limit argument, we can then extend the Wick product to more general random 
variables (see [T7] for further details). In this paper, we will take the limits in the L 2 (Q) 
topology. 

For / G H we define its exponential vector S(f) by 

S(f) := e<*C/> = exp U(f) - MMi\ = eX p U(f) - (h(f)f 



OO oo 1 

z — ' n\ z — 4 n\ 

n=0 n=0 



In a similar way we can define the complex exponential vector of / by (z denotes here the 
imaginary unity) 

e-M/) = exp L W) + Wi\ = l - ]h {fr. (87) 

With these notations in hand, an important property of Wick product is the following 
relation: for any two elements / and g in Ti, we have 

£{f)*£(ff)=£(f + g), (88) 
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an analogous property for the complex exponential vector being also satisfied. 
We now state a result which is a generalization of (TT| Proposition 4.8]. 

Proposition 6.1. Let F £ DomD k and g £ H m . Then 

(1) F o Ik{g) is well defined in L 2 (Q). 

(2) Fg £ Dom5 ok . 

(3) FoI k (g) = 5° k (Fg). 

Proof. Let F £ DomD fe . This implies that F admits the chaos decomposition F 

E^=0 J n(/«)' With 



J> fc n! ||/ B | 



< oo. 



(89) 



n=0 



Define then Fn = ^2 n= o In(fn)- Consider also g £ T-i® k . In order to check (1), we shall see 
that the limit in L 2 (Vt) of Fj^oI k (g) exists, as iV — > oo. But Fj^oI k (g) = X^ =0 In+k{fn®g), 
and the limit in L 2 (Q) of this expression exists if and only if 



J2(n + k)\\\f n ®g\ 



< oo. 



n=0 



This last condition is clearly satisfied thanks to (189]) . 

Now we will prove our claims (2) and (3) for F with a finite chaos decomposition and 
g = 9\ ® ■ ■ ■ ® gki with gi £ TL, for i = 1, . . . , k. More precisely, we will see that for any 
G £ S the following relationship holds: 



E 



{FoI k {g))G =E[(Fg,D k Gy, 



(90) 



This will be done by an induction argument. For k = 1, this is a consequence of |17[ 
Proposition 4.8] since Fo I\{g) = S°(Fg). Suppose now that (19"0|) is satisfied for k = K. 
Therefore, 



E 



[FoI K+l {g)]G 



E 
E 



(F 0/1(01) 04(52 
{(FoI 1 {g 1 ))g 2 ®- 



■®9 K+ i))G 

g K +i > D K G) n ®K 



where in the last equality we have used that F o Ii (gi ) has a finite chaos expansion and 
the induction hypothesis. The last expression can be rewritten as 



E 



(F o h(gi)) D^... mK+ G 



Since D^...^ gK+i G £ S, we can apply the case k — 1 to the above expression and we 
obtain 



E 

=E 



[F«I K+1 (g)]G 



E 



{Fgi , D 



K 

9 2 i 



G) 



H 



FD 



K+l 



G 



E 



(Fgi 



g 



E 



D K+1 G) 



FD 1 (D 

91 V 



K 

g 2 i 



K+l ' 



which finishes our induction procedure. Thus, fl90j) is satisfied for F with a finite chaos 
expansion and g a tensor product of elements of "H. 

To extend the result to a general F £ Dom(D /c ) and g £ "H® fc , we first consider the 
case F £ ~Dom.(D k ) and g = g\ (g) • • • ® g k . In this situation, identity (190]) is a consequence 
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of the fact that this relationship holds for Fn = Yln=o^n{fn) defined above and of the 
part (1) of the proposition. Finally, for a general g G using the fact that both sides 

of f )90|) are linear in g, we can generalize this identity to g belonging to the linear span 
of elements of the form gi <£>■■■ <S> gk, which is a dense subspace of "H® fc . So if g 6 H® h 
and {g M }MeN is a sequence of elements of this linear span of tensor products such that 
g M — y g in H® k , one can easily see (by using (I89p ) that 

FoI k (g M )^FoI k (g), 

as M — y oo in L 2 (Q). Since 



E 



{FoI k {g M ))G = E[(Fg M , D k G) 



the proof is completed by a limiting argument. □ 

6.2. One-dimensional case. In order to simplify a little our presentation, we first show 
the identification of 5°(Vf(x)) with lim|n 3t |->.o S UsU<> when d — 1, that is when x is a one- 
dimensional Gaussian process satisfying Hypothesis 15.11 The first step in this direction is 
a general formula for Wick products of the form G(X) oY op , where X and Y are elements 
of the first chaos (see Section ISTTl for a definition). Notice that the proof of this proposition 
is deferred to the Appendix for sake of clarity. 

Proposition 6.2. Let g,h G H and let G : R — > R be differentiable up to order p such 
that all its derivatives G^' are elements of L r {jj, g ) for any j = 0, . . . ,p and for some 
r > 2, with fig = jV(0, ||y||^). Define X = h(g) and Y = I\{h). Then the Wick product 
G(X) o y op can be expressed in terms of ordinary products as 

G(X) o Y op = G{X)Y P 

( 1 \m+l n \ 

+ £ 2-mHv - 2m- fV ^ W [E(XFtf ^1^' " (91) 

0<l+2m<p >' 

Example 6.3. In order to illustrate the kind of correction terms we obtain, let us write 
formula f }9T|) for p = 1,2, 3: 

G{X)oY = G{X)Y - G'(X)E(XY) 

G(X) o F° 2 = G(X)Y 2 - G(X)E(Y 2 ) - 2G'{X)E(XY)Y + G"{X) [E(XY)] 2 

G(X)o7° 3 = G{X)Y 3 -3G{X)E(Y 2 )Y + 3G'{X)E{XY)E{Y 2 ) 

-3G'(X)E(XY)Y 2 + 3G"(X) [E(XY)] 2 Y - G"\X) [E(XY)f . 

We are now ready to state our representation of the Skorohod integral by Riemann-Wick 
sums: 

Theorem 6.4. Let x be a 1-dimensional centered Gaussian process with continuous co- 
variance function fulfilling Hypotheses I5.il and \5.2\ and assume that x also satisfies Hy- 
potheses \L1\ Let f be a function in C 27V (R) such that fW verifies the growth condition 
(GC) for k = 1, . . . ,2N. Then, the Skorohod integral 5°(Vf(x)) (whose existence is en- 
sured by Theorem \5.13\) can be represented as a.s. — limn s t-i.o S UsU<> ; where S Ust, ° is defined 
by 

n-l N 1 
i=0 k=l 
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Moreover, we have 



**(/'(*)) = / f'Mdx p 



f(x p )R'dp 



(92) 



Proof. As mentioned at the beginning of the section, our main task is to compute S 11 ^' 
in terms of ordinary products. This will be achieved by applying Proposition 16.21 to each 
term in the above sum, with G — f , X — x ti and Y = x£. t . 1 = x ti+1 — x ti . 

To this end, notice first that the integrability conditions on / required at Proposition ^. 21 
are fulfilled as soon as condition (GC) (see Definition l5.1ip is met. Fix then i G {0, . . . , n — 
1}, recall that we set X = x^ and Y = x^. t . l? and consider the quantity S l k := ^/^(xtjo 



\o,k 



titi + l I 

N 



A direct application of Proposition 16.21 yields 

( 1 \l+m 



E S * = E E 



2 m m\l\(k - 2m - /)! 



/(*+*) (X) [E(ir)f [E(F 2 )] m F 



k—2rn—l 



k=l k=X l+2m<k 

Making a substitution q = k + I and I + m — u, this expression can be simplified into 

/+m<?/2 

m ■ ■ r E(y 2 ) 



E s i 



fe=l 



2 m m\l\(q - 21 - 2m)! 



3=1 

27V 

E^ (9) ( X )E (, 

<?=1 «=0 v ^ 7 l+m=u 

27V [9/2] 

E 



•1)" 



(g-2u)!u! 



E(jy) + 



E(F 2 



-2« 



Moreover, recalling again that X = x ti and Y = x ti+1 — x ti , it is easily seen that 

E(Y 2 ) 1 



E(XY) + 



e« +i )-e(4; 



(93) 



(94) 



where we recall that 5R ti t i+1 stands for Rt i+1 — -R< r Therefore, summing now over i G 
{i, . . . , n — 1} we get 

n-l AT 27V [q/2] 



-1" 



n-l 



s*w j^^s* EE( g _ 

g=l U =0 V ^ ' 1=0 



(95) 



i=0 k=l 

where the quantity 7^' u is defined by 

We now separate the study into different cases. 

Case 1: If u = 1 and g — 2u = (namely q — 2), then 

n-l _. n-l -. pt n—l 



(96) 



i=0 



i=0 i=0 

where in the last equality we resort to the fact that R p is absolutely continuous (see 
Hypothesis 15. ip . From this expression, by a dominated convergence argument one easily 



\— \n— 1 rrQ.U 

gets hm^oo 2^ - =0 77' 



Ulnx^dp. 
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Case 2: If u > 2 or u = 1, q — 2u > 1, then lim^oo Y^=o = 0- Indeed, recalling 



-1-1,11 

n—l „f n—l 



definition (196]) of Tf' , we observe that 



IL _L IL _L 

i=0 * ™ ^ s i=0 

In the right hand side of the above inequality, it is now easily seen that 

hm max { |x^ +i I"" 2 ", I^^J"" 1 } = 0, 

n->ooO<J<n-l L + J 

while the integral term remains bounded by C f* \R! \ dp, which is bounded by assumption. 
This completes the proof of our claim. 

Case 3: If u = 0, then 

n-1 n-1 JV n-1 27V 

E^* = EE >'(*.,) «..,«)" + E E >(*<.)«,«,+,)'• 

i=0 i=0 9=1 y ' i=0 g=/V+l V ' 

Thus Theorem 14.51 asserts that the first sum above converges to f f'(x u )dx u , while it is 
easy to see that the second sum converges to 0, thanks to the regularity properties of x. 

Plugging now the study of our 3 cases into equation (195]) . the proof of our theorem is 
easily completed. 

□ 

6.3. Relationship with existing results. Several results exist on the convergence of 
Riemann-Wick sums, among which emerges |27| . dealing with a situation which is similar 
to ours in the case of a one-dimensional process. 

In order to compare our results with those of [27], let us specialize our situation to the 
case of a dyadic partition of an interval [s, t] with s < t (the case of a general partition is 
handled in [27], but this restriction will be more convenient for our purposes). Namely, 
for n > 1, we consider the partition ITJ = < k < 2 n }, where t% = s + k(t — s)/2 n . 
For notational sake, we often write tk instead of t^. We shall also restrict our study to 
the case of a fBm B, though [27] deals with a rather general Gaussian process. 

Let us first quote some results about weighted sums taken from [T2| HP] : 

Proposition 6.5. Let B be a one- dimensional fractional Brownian motion, whose covari- 
ance function is defined by Let g be aC 4 function satisfying Hypothesis (GC) together 
with all its derivatives. 

(i) For n > 1, set 



2' 1 — 1 
^ 1 



~2nH 



v^ 2) (g) = Y,9(B tk : 

Then if 1/4 < H < 3/4, we have 

C - hm ^"^(j) = a H f 3 g{B s ) dW s , (97) 

where £ — lim stands for a convergence in law, an is a positive constant depending only 
on H, and W is a Brownian motion independent of B. 
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(ii) For n > 1, set 

2 n -l 

K (3) (s) = 5>(^J (b* 

k=0 

Then if H < 1/2 we have 

L 2 (n) - lim n' H - l V^\g) = -\ f S g'(B s )ds. 

We can now recall the main result of [27J, to which we would like to compare our own 
computations, is the following: 

Proposition 6.6. Let B be a one- dimensional f Bra with Hurst parameter 1/4 < H < 1/2 
and f be a C A function satisfying Hypothesis (GC) together with all its derivatives. For 
< s < t < T , consider the set of dyadic partitions {n™ t ; n > 1} and set 

S^ = ^f(B tk )oBl tk+i . (98) 

fc=0 

Then S n,<> converges in L 2 (Q) to 5°(f'(B)) (which is the Skorohod integral introduced at 
Theorem \5.13\) . 



Proof. Our aim here is not to reproduce the proof contained in [27], but to give a version 
compatible with our formalism. We shall focus on the case 1/4 < H < 1/3, the other one 
being easier. 

Consider first < u < v < T. According to Example I6.3[ we have 
f'(B s ) oBl = f\B u ) Bi„ - ^f"(B u ) E [B u B*J 

= f'(Bu) Bl - \f"{B u ) [v 2H - u 2H ] + \f'{B u ) \t - s\ 2H . 
In a rather artificial way, we shall recast this identity into 

3 1 

f(B s ) o Bl = -/ 3 \B U ) (Bl) 1 - Rl - Rl, (99) 

with 

Rl = \f"{B u ) [(BL) 2 -\v- u\ 2H ] , and R 2 UV = \f^{B u ) (B^) 3 



1 

2 

Plugging f l9"9~j) into the definition of S n,<> , we thus obtain 

2 n -l 



S">* = - \ £ f\B tk ) [til, - 4 H } + \v?\P) + Jk (3) (/ (3) ). (ioo) 

k=0 



Now, invoking Proposition 16.51 it is readily checked that both Vn 2 \f") and Vn(f^) 
converge to in L 2 (Q) as n — > oo. Hence 

L 2 (n) - lim = J st (f'(B) dB) -H f f"(B u ) u 2 "" 1 du, 
which ends the proof. 

□ 
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The aim of the computations above was to prove that the results of |27] do not contradict 
ours for H > 1/4. Note however the following points: 

(i) Having a look at Proposition 16.61 one might think that the first order Riemann- 
Wick sums S n,<> are always convergent in L 2 (Q). However, when H < 1/4, relation (j97p 

(2) 

still holds true. This means that the term Vn (/") appearing in equation f llOOp is now 
divergent, due to the fact that 2H — 1/2 < 0. The same kind of arguments also yield the 
divergence of ^^(f^) in (UOOp . It is thus reasonable to think that first order Riemann- 
Wick sums will be divergent for H < 1/4, justifying our higher order expansions. 

(ii) In light of Proposition 16. 6[ it is however possible that expansions of lower order than 
ours are sufficient to guarantee the convergence of sums like S Uat,<> in Theorem 16.41 We 
haven't followed this line of investigation for sake of conciseness, but let us stress the fact 
that almost sure convergences of our Wick-Riemann sums are obtained in Theorem 16.41 
and Theorem 16.81 (for any sequence of partitions whose mesh tends to 0), while only L 2 (Q) 
convergences are considered in Proposition 16.61 

(in) It is also worth reminding that we aim at considering a general rf-dimensional Gauss- 
ian process, while [T2| [27] focus on 1-dimensional situations. It is an open question for 
us to know if the methods of the aforementioned papers could be easily adapted to a 
multidimensional process. 

6.4. Multidimensional case. We shall now give the representation theorem for Skoro- 
hod's integral in the multidimensional case. Technically speaking, this will be an elabo- 
ration of the one-dimensional case, relying on tensorization and cumbersome notations. 

We first need an analogous of Proposition 16. 21 in the multidimensional case, whose proof 
is also postponed to the Appendix. To this aim, let us introduce some additional notation: 
given gi, . . . , g d EH define g — (gi, . . . , g d ) and denote by \ig the law in M. d of the random 
vector (ii(s-i), . . . , h{g d )). 

Proposition 6.7. Using the notations introduced above, let g±, . . . , g d , hi..., ha G "H. 
Consider the random variables X\ = ii(<7i), . . . , X d = Ii(gd), and Yi = I\(hi), . . . , Y d = 
Ii(hd)- Suppose that Y\, ■ ■ ■ ,Y d are independent and also that Xj and Yk are independent 
for k 7^ j. Let p = (pi, ■ ■ ■ ,Pa) be a multiindex and set \p\ = Y^j=iPj- Assume that 
G e C^(R d ) is such that d a G G L r (fi- g ) for any multiindex a = (cki . . . , a^) and for some 
r > 2, with 0!^ < p/., k — 1, . . . , d. Then G(Xi, • ■ ■ , X d ) o Y° P1 o ■ • • o Y^ Pd is well defined 
in L 2 {VL) and the following formula holds: 

G(X h --- ,X d )oY^o---oY; pd = ••■ & u - M G{Xi,--- ,X d ) 

x n (nXtYt)),k {nY ' )r yr " 2m *"'*] ■ (101) 

where d^ 1, "'^ d denotes dJ )^ . 

As in the one-dimensional case, the proposition above is the key ingredient in order to 
establish the following representation formula for Skorohod's integral: 
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Theorem 6.8. Let x be a d- dimensional centered Gaussian process with continuous co- 
variance function fulfilling Hypotheses 15.11 and \5.2\ and assume that x also satisfies Hy- 
pothesis \l.l[ Let f be a function in C 2N (R d ) such that d a f verifies the growth condition 
(GC) for any multiindex a such that \a\ < 2N . Then the Skorohod integral 5°(V/(x)) 
(whose existence is ensured by Theorem \5.13\) can be represented as a.s. — limn st ^o S nst ' , 
where S Ust,<> is defined by 

n-l N 

5 rw = J2 J2 ^ „/Xo ; ) o *lu + , o • • • o (4). 

i=0 k=l 

Proof. We mimic here the proof of Theorem 16.41 decompose first S 11 "*' into Y^=o Z~2u=o 
, , „■ , where 

^-^H k7d=« 31,—,Jd J 

SL.,u = JT^^JMV* ■ ■ ■ >*u(d)) o (x^l))* 1 o ■ ■ ■ o (xj, +1 (d))^ • 

jl- ' ' Ja- 

For a fixed i G {0, . . . , n — 1} and fc G {0, . . . , d}, set now 

A fc = x u (k), and Y" fc = xJ t . +i (A;). 

As in the one dimensional case, it is readily checked that if d a f satisfies the growth 
condition (GC) for any \a\ < 2N, then the integrability conditions of Proposition 16.71 are 

also fulfilled for x ti = (x^(l), . . . ,x ti (d)) = (^i(lpti))' ■ ■ ■ ' -Ml[o *<)))■ This anows to write 



JV 

Y Y s 'i .i.; Y Y Y "' Y d h+ju-,i d +uf( x i,--- ,x d ) 

u=l jiH hj'd=« u=l iiH hjd=« ii+2mi<j'i ld+2m d <j d 

d ( -I \(w,k+lh) 

TT XJJ. nx k Y k )) h (E(Y k 2 )) mk Yi k ~ 2mk ~ h 



Making substitution + j k = q k for k — 1, 2, • ■ ■ , cZ, or jj. = q k — l k , the condition 
Zfe + 2mfc < jfc can be written as 4 + m k = u k with < u k < q k /2 and therefore the same 
kind of manipulations as in f l§5j) yield 



N 

Y Y s L,u - Y d qu ..., q j(x 1 ,..-,x d ) 

u=l jx+-+j d =u l<qi+---+q d <2N 



X 

A- 
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Furthermore, like in the proof of Theorem 16 A\ we have E(Xfclfc) H = 5R tl t l+1 /2. 

Summing over i G {i, . . . , n — 1} we thus end up with 

n—1 d 

^ = e e ^••^/(^( i )r--,^(d))n^(xi^ 1 (A ; )) ft 

t=0 l<gi+-+9 d <JV fc=l ^* fc '" 



n-1 



i=0 jV+l<gi+-+g d <2JV fc=l 
1 n—1 (i 

"2 EE^V(^(1), • • ■ - Ru) + 

i=0 k=l 

d 

:= 0?" + + E ©3? + @ nst • 
fc=i 

In the last sum, the variables correspond to the 3 cases we have distinguished in 
the proof of Theorem 16.41 G^ st denotes the sums in which the Uk are equal to and 

1 < qi H \-qd < N; 0" st are the terms with u k = and A" + 1 < <?i H V qd < 2iV; 

3 ^ corresponds to the terms with q^ = 2, qj = for all j ^ k and — 1 (so that Uj = 

if j j£ k). Finally, Q Ust denotes the sums with either u\ + ■ ■ ■ + u<i > 2 or some Uk — 1 
(and = for j ^ k) but > 3. Referring again to the proof of Theorem 16 A\ it is then 
easy to argue that nst and 0^ st converge to 0, 0f s * converges to / s *(V/(x p ), dx p )^d, and 
Efc=i Q 3,k converges to -\ J* Af(x p )R' p dp. 

□ 

7. Appendix 

In this Appendix, we prove Propositions 16.21 and 16. 71 For this, we will need the following 
analytical lemma. 

Lemma 7.1. Let fi be a finite measure on (M d , B(M. d )) and let G G C p (M. d ) be such that 
d a G G L r (fi) for some r > 1 and any multiindex a such that \a\ := Ylj=i a j — P- Then, 
there exists a sequence (G n ) n eN such that 

(1) Each G n is a trigonometric polynomial of several variables, that is G n (xi, . . . ,Xd) = 
Sfinite* 3 ™! h xi+ ' l d Xd , and where ld and are real numbers. 

(2) We have \im. n ^ f00 d a G n = d a G in L r (p) for any a such that \a\ < p. 

Proof. This lemma is folklore, but we haven't been able to find it in any standard text 
book. For this reason and for the sake of completeness, we give here the main ideas of its 
proof. First, given G G C p {R d ), there exists a sequence of C°°(M a! ) functions with compact 
support that converge, jointly with their derivatives, to G in L r (fi). So, one only needs 
to approximate a function G G C°°(K d ) with support contained in a rectangle of M. d , say 
K. Moreover, given e > 0, we can suppose that n(K c ) < e. For a such function, consider 
its Fourier partial sums on the rectangle K that converge uniformly, jointly with their 
derivatives to G and its derivatives. Since these partial sums are periodic functions with 
the same period, their sup-norm on all M. d is the same that the sup-norm on the compact 
K. With these ingredients, the result is easily obtained. 
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□ 

Proof of Proposition ^. 2l We start with G(x) = e l ^ x , for an arbitrary (GK, which means 
that we wish to evaluate the Wick product e^ x o Y op . 

Recall that X = I\{g) and Y = Ii(h). For £,77 G R, consider the random variable 
M(£, V ) = exp f-^E(X 2 )^) S (z& + 1*7/1) = exp ( -^E(X 2 )) £ (1&7) o £ (177/1) 



2 



exp o £ (177/1) = — T ex P « X ) F ° P 



p=0 ^ 

where we have invoked relation (I88|) for the second equality and relation (1ST)) for the last 
one. It is thus obvious that e l ^ x o Y op can be expressed as 

p\ 

— x the coefficient of rf in the expansion of M(£, 77). 
We now proceed to this expansion: we have 

M(fij) = ex P | ! ^x + ! i ! y + ^E(y 2 ) + ^E(xy) 

- ^E^E^tWE^TOf. 

fc=0 m=0 Z=0 



1- ^ \E(XY)] 1 \ E (Y 2 )] m YP- 2m - 1 



-2m-/ 



Hence, by computing the coefficient of ?7 P in the above expression, it is easily checked that 

^-2m-l^lp\ 

2 m m\ l\ (p-2m-l)\ 

<* Y p + V ( —e*A (~ 1 ) m+f P ! \E(XY)} 1 \E(Y 2 )] m Y p ~ 

+ 2^ \ dx i e )2™m\l\(p-2m-l)\ [{ )l[{ )1 

which is the desired formula fl9Tj) for G(a;) = e^ x . 

Let us now see how to extend this relation to a more general function G. By linearity, 
we first obtain the result for any trigonometric polynomial G. Now, let G be such that 
qU) ^ L r ([ig) for any j — 0, . . . ,p and some r > 2. By Lemma [77T| there exists a sequence 
(G n ) nS N of trigonometric polynomials such that 

G® -)• G7 (i) in L r (/i 9 ) for any j = 0, . . . , p. 
This implies that G(X) G Dom(D p ) and that 

D>G{X) = G {j) (X)g® j for any j = 0, . . . , p. 
Indeed, G n (X) G S and D^G n (X) = G < i ) {X)g® j for any j = 0, . . . ,p. Moreover, since 

E[\G<p(X) - G^(X)\ r ] = \\Gtf - G®f LrM 

we have that 
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Using that the D^' are closed operators we obtain that G(X) E Dom(D p ) and that 
D j G(X) = G^{X)g® j for any j — 0, . . . ,p. In particular, owing to Proposition 16 . 1 1 we 
have that 

G{X) o Y op = G(X) o I p {h® p ) = 5 op (G(X)h® p ). (102) 

Let us go back now to our approximating sequence {G n ) n ^. It is readily checked that 
relation fl 1 2 [) also holds for any G n . Thus, putting together the relation G n (X)oI p (h® p ) = 
8 <>p {G n {X)h® p ) with equation (}9Tj) for a trigonometric polynomial, we get that 

5 0p (G n (X)g® p ) 

= G n (x)Y p + £ [ J _ 2 P m _ v Gg W mxY)f [E(y 2 )]-y^(i03) 

Since -> GW(X) in L r (fi) with r > 2 and the Y k - 2m ~ l belong to all the 

the right-hand side of f ll03p converges in L 2 (f2), as n — > oo, to 

G{X)YP+ £ , j _ ^ _ (X) [E(IF)]' [E^f . 

0<Z+2m<p V ^ ' 

Finally we obtain the general case of equation (l9Tj) by taking limits in both sides of 
equation fll03[) and by resorting to the closeness of the operator 5 op . 

□ 



As in the previous section, the extension of Proposition 16.21 to the multidimensional 
case is now an elaboration of the previous computations relying on some notational tech- 
nicalities. 

Proof of Proposition 1 6. 71 As for Proposition 16.21 we first consider G(x) = e l ^' x \ where 
x = (x\, . . . , Xd) and £ = (£i, • • • , £d) are arbitrary vectors in R d . The extension of the 
formula to a G satisfying the general integrability conditions of our hypotheses is then 
obtained following the same approximation scheme as in the one-dimensional case, and is 
left to the reader for sake of conciseness. 

In order to treat the case of G(x) = e l ^' x \ set 

/ d 
M(£,7?) = exp (i<£, X)) oexp i( V , Y) + - ^%F fc 2 

V k=l 

Along the same lines as for Proposition 16.21 one can then identify e*^' x ^ o Y° P1 o • • • o Y^ Pd 
with p^+" p +p d x the coefficient of t] Pl ■ ■ ■ rf d d in the expansion of M(£, rj). Moreover, thanks 
to relation (|88j) and invoking the fact that Xj and Y k are independent for k ^ j, we get 
that 

i X ZkX k + z J] + - vMYk) + X ZkVk~E(X k Y k ) 

k=l k=l k=l k=l 
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Expanding now the exponential according to formula (|87j) . we end up with 

d t (p k -2m k -2l k ) 



M(t,v)= E E •• E n{(^) /fce< 

Pl,-,Pd=0 h+2m 1 <p 1 ld+2m d <p d k=l 



e.X k 



2 m *m k \l k \(p k - 2m k - l k )\ 



x 



(E(x k Y k )t (E(n 2 )) mfc y, 



2\\ m k vPk -2m k -l k 
k 



Pi 

m 



Vd 



Taking into account the fact that Y[t=i(' l £k) lkel ^-' k=1 ^ kXk — <9' 1 >---> Zd e l ^= 1 ^ cXfc , our for- 
mula fllOip is now easily deduced, which ends the proof. 

□ 
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